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ABSTRACT 



Numerical integration and Monte Carlo techniques are 
used in the development of several models in order to 
determine the effect on probability of random detection 
of a merchant ship using speeds up to 90 knots by a 10 
knot submarine patrolling a back-and-f orth barrier. A 
definite range law for detection is assumed. Individual 
encounter models are developed for ship tracks which cross 
the midpoint of the submarine patrol line at various angles. 
The models are extended to include the assumption of a 
normal-distribution of crossing points. Computer programs 
of the models, written in the FORTRAN IV language, are in- 
cluded. The results are applied in a numerical example. 

It is concluded that while increases in ship speeds do 
result in a substantial decrease in probability of detection 
by a submarine in the case of a single barrier transit, a 
speed advantage alone when applied to a typical transit of 
the North Atlantic will not appreciably decrease the overall 
detection probability. 
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I. 



INTRODUCTION 



A. BACKGROUND 

1 . High-Speed Ship Development 

New principles and concepts in ship hull design 
have made possible the development of ships capable of 
speeds heretofore considered impossible. In particular, 
the advent of surface-effect ship (air cushion and captured 
bubble), hydrofoil, and catamaran hull designs, coupled 
with advances in nuclear power, gas turbine, and water- jet 
propulsion systems, will have revolutionary effects on 
naval strategy and tactics, and world shipping patterns. 
That serious consideration to high-speed merchantmen is 
being made is evidenced by the fact that the U.S. Navy and 
the Maritime Administration share funding in the Joint 
Surface-Effect Ship Program which was established in 1966 
to determine the feasibility of 4,000 to 5,000-ton surface- 
effect ships capable of speeds in excess of 80 knots. Bell 
Aerosystems, Inc., and Aerojet-General Corp. have each 
received $1.5 million as a start on a 3-year program to 
build an experimental 100-ton surface-effect ship [1] . The 
fact that high-speed merchant ships which can enjoy large 
speed advantages over submarines are indeed in the offing 
indicates a need for reexamination of the abilities of 
future submarines to detect and kill these ships. This 
thesis makes a start by examining the effect of ship- 
submarine speed ratios on detection probability 
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resulting from random encounters by submarines on fixed 
patrols . 

a. Literature Search 

As might be expected, the literature dealing 
with high-speed ships is relatively sparse; those papers 
that do make reference to ships of the future concern 
themselves with the hydrofoil small-ship types [2], or with 
the so-called super merchant ships which have such features 
as large (500,000-ton) displacements, and underwater 
tankers of speeds up to 40 knots [3]. Various Navy studies 
in the past which have attended to the question of ship 
speeds versus vulnerability have been limited in scope to 
surface ship speed ranges up to only 35 knots [4, 5]. One 
study in particular concluded that "...When ships are sailed 
independently. . . (and) when cost is considered there is only 
a marginal gain in going from 20 to 24 knots...." /!/ . The 
cost consideration at that time (1955) involved studies of 
conventional hulls. An earlier study [6] showed a decrease 
in percentage of independent merchant ships sunk from 10-12 
per cent when ship speeds were approximately equal to 
submerged submarine speeds (4-8 knots) to 1-2 per cent when 
ship speeds were double that of the U-boat submerged speeds. 
When hull-mounted sonars on merchant ships were considered 
[7], maximum optimal sustained speeds for independent 
vessels were in the range of 20 to 24 knots. 

In summary, neither the search of available 
studies [8], operations research papers, nor computer 
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searches for documentation [9] at the Defense Documentation 
Center, Washington, D.C., and at the U.S. Naval Postgraduate 
School, Monterey, California have revealed any literature 
regarding high-speed ships, and in particular the effect of 
merchant ship speed increases beyond 35 knots on vulner- 
ability to submarines. Indeed, all facets associated with 
these high-speed ships should be examined, e.g., noise 
radiation in water, the submarine fire control problem, 
effect of torpedoes, etcetera. 

2 . Vulnerability of Convoys and Independent Ships 

As was clearly evidenced during World Wars I and II, 
convoys and independent merchants were quite susceptible to 
attack by submarines. Allied shipping losses from the 
German U-boat campaign during the first five years of World 
War II provide sufficient evidence of the extreme vulner- 
ability of merchant shipping [10, 11]. Considering the 
improvements in both submarine and submarine weapon system 
capabilities since this era, the shipping losses previously 
experienced can be regarded as the minimum losses that may 
be expected in a future war [10]. Furthermore, it is a 
foregone conclusion that in the opening stages of a future 
general war escort numbers will be insufficient to provide 
adequate convoy protection for all merchant shipping [12, 
13]. It can be assumed, then, that the initial phases of 
a future war will see an aggressive enemy submarine force 
presented with many unescorted targets which can be attacked 
with impunity [10]. 
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To counter the possibility of sustaining heavy 
merchant losses at the outset of a future war, passive 
defense measures such as installation of sonar on merchant 
ships, routing to minimize losses, convoying, and increas- 
ing merchant ship speeds have been considered [12]. This 
thesis considers the passive defense measure of increased 
merchant ship speed. 

3 . Future Capabilities of Submarines 

The Type XXI submarine produced by the Germans in 
1944 influenced basic submarine design until the appearance 
of the nuclear-powered submarine in 1954 [14]. Since the 
introduction of the nuclear-powered submarine, great tech- 
nological strides have been realized in the areas of sub- 
marine hull construction, propulsion, sonar, torpedo 
weapons systems, and fire control systems. These advances 
have contributed significantly to the threat to merchant 
shipping posed by submarine forces of potential enemies. 

The modern nuclear-powered submarine can equal or exceed 
the speed of its surface-ship target or pursuer, and can 
remain submerged for days. Together with this endurance 
capability, the nuclear-powered submarine is able to main- 
tain a high sustained submerged speed in a relatively 
quiet environment. "...The limiting lines of approach to 
convoys and task forces steadily widen with increasing 
submarine speed until there is no sector safe from attack 
by a 30-knot submarine...." /2/ This quote was taken from 
a publication dated April 30, 1962; since this time a great 
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deal of research and development in the area of surface- 
effect ships has taken place. The Joint Surface-Effect 
Ship Program was initiated to consider large surface- 
effect ships capable of high speeds. It is conceivable 
that in the forseeable future surface ships with speeds 
far exceeding the capabilities of the nuclear-powered 
submarine will take their place in the surface ship com- 
munity. Submarine speeds are not expected to exceed their 
present maximum for some time to come , thus permitting 
the high-speed merchant ship to enjoy a comfortable ad- 
vantage in speed. 

4 . Safety of Independents 

Heretofore convoys have been recommended as the 
best alternative for moving merchant vessels across the 
oceans when a submarine threat exists [15]. This is due 
primarily to the fact that for slow moving (less than 15 
knots) merchants, safety lies in numbers. The measure of 
effectiveness for analysis of merchant shipping has been 
percentage lost of the total number of ships in the 
convoy [16]. Studies completed during World War II indi- 
cated that approximately the same number of ships were lost, 
per attack, from large convoys as from small convoys. Hence 
the conclusion was to not only continue convoying, but to 
increase the size of the convoys as well. 

With the advent of high-speed surface-effect ships, 
the problem of independent sailing versus convoying should 
again be considered. In a study conducted by the 
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Antisubmarine Operations Research Group Tenth Fleet in 
December 1943, the advantages and disadvantages that might 
result from independent sailings were discussed: "....Fast 

merchant vessels, with the exception of those carrying high 
priority cargo, have been able to continue sailing inde- 
pendently along the coast (Atlantic) much of the time." /3/ 
Here is a direct reference to a merchant ship speed advan- 
tage over the submarine. 

Other advantages of independent sailings are of 
two kinds: first, the independent merchant can operate more 

efficiently and hence deliver more cargo in a given time; 
second, the reduced requirements for escorts releases these 
surface, air, and submerged craft for offensive roles [11]. 

5 . Probabilities of Encounter 

The problem of random encounter between two units, 
sometimes referred to as the searcher and the target [17, 

18, 19, 20], is examined in most of the classical operations 
research literature. Operations Evaluation Group reports 
and studies (which constitute a large part of the so-called 
"classical literature"), and more recently by search theor- 
ists and game theorists. 

Probabilities of encounters, or detections resulting 
from applications of search plans , are discussed in much of 
the literature, ranging from Kimball and Morse [16] and 
Koopman [18], through the summary reports of the Operations 
Evaluation Group [14] and applications studies [2] to the 
modern search theorists, such as Pollock [19, 20]. 
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Data concerning encounters versus ship position 
displacement from reported submarine contacts [21] and 
frequency of encounters versus submarine position predic- 
tions [22] have been gathered and analyzed during the 
Second World War. Applications of game theory for the 
selection of convoy routes which would minimize ship- 
submarine encounters were proposed in a study by Danskin 
of the Operations Evaluation Group in 1954 [15], which 
later appeared as a paper in Operations Research [23], The 
Random Encounter with no involved search effort has been 
studied by several authors; Koopman [17, 18] considered the 
random encounter of an observer progressing at a constant 
velocity within a uniform distribution of targets of con- 
stant speed, as well as within a circular-normal distribu- 
tion of targets. The problem of random encounter between 
a surface force progressing at constant velocity with 
another unit travelling at a fixed speed, the position and 
heading of which are independently and normally distributed, 
has been examined by Dobbie in 1945 [24]. 

Literature concerning the effect of speed on the 
random encounter is limited almost exclusively to government 
studies and reports. Kittel, in 1944, analyzed data con- 
cerning sinkings and sightings by submarines to give an 
early indication of the effect of speed on the safety of 
ships [6]. The results were discussed by Sternhell and 
Thorndike in 1946 [25], and by Winston in 1955 [4], These 
results, in turn, were summarized by the Operations 
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Evaluation Group in 1957 [12]. More recently, the general 
problem of the effect of speed on the vulnerability of 
independent merchant vessels was examined by Neufer in 
1961 [4]. Indeed, Neufer appears to have been the first 

to attempt to determine the explicit survival probability 
of a merchant vessel in transit from one port to another as 
a function of its speed. 



a barrier patrol has been examined by at least two authors. 
Koopman , in his Search and Screening, examined the problem 
using the underlying assumption of a uniform distribution 
of points of intersection of the target's track with the 
back-and-forth barrier patrol line. /4/ In particular, 
Koopman showed that, with a definite range law, the 
probability of detection by the back-and-forth barrier 
patrol is given by 



Finally, the probability of encounter by a unit on 




/ 2 ’ 2 

[A - (/l/6 + 1 - l)/2] /A (A + 1) 



for 6> 1/ (2/A (A + 1)' ) 



\ 1, for 6 < 1/ (2/A (A + 1)' ) 



where 



A = ratio of the barrier patrol line length to 
sweep width 



<5 = ratio of target speed to barrier unit patrol 
speed 
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Further, Koopman showed that, for all values of A, the 
back-and-forth barrier patrol is preferred to the crossover 
("figure eight") type patrol for all values of 6 l 1. /5/ 
Van Train [26] used a war gaming method in a limited way to 
determine the probability of detection of a transitor 
intercepting a multi-unit, single or multi-line fixed 
barrier, in which each submarine was constrained 

1) to remain in the center of an assigned zone, and 

2) to conduct intercepts using headings parallel to the 
barrier line. The submarine's course, and position on a 
"starting line" were assumed to have been normally- 
distributed . 

B. THE PROBLEM 

The basic encounter problem is that of determining the 
probability of detection of a surface ship by a submarine 
on an assigned back-and-forth barrier patrol, such as that 
discussed by Koopman [18] . The following assumptions were 
made for the models: 

1 . Assumptions for the Basic Random Encounter Model 

a. The submarine patrols a geographically-fixed 
back-and-forth barrier centered on a known surface ship 
track or convoy route . 

b. The submarine patrols at a constant speed, v. 

c. A definite range law for detection is used. 

d. The submarine patrol length, L, is four times 
the detection radius. 
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e. The submarine position at any time, x, is from 
a uniform distribution. 

f. The point of intersection of the ship's track 
and the patrol line (or extension of the patrol line) is 
from a normal distribution, with standard deviation of 1/6 ^ 
of the length of the. submarine patrol line. 

C. METHODS OF SOLUTION 

Two methods of solution to the models were devised; the 
first was a numerical integration technique, i.e., simula- 
tion by systematic sampling. The submarine's initial 
position was successively incremented by 500 yards (0.25 
nautical mile) and the minimum lateral range (CPA) between 
the submarine and the surface vessel was computed by 
numerical vector analysis. The basic model considered 960 
initial submarine positions for each speed ratio considered; 
twelve ship-submarine speed ratios , ranging from 1 to 9 were 
examined. The basic models also examined nine angles of 
intersection of tracks, ranging from 90 down to 30 degrees. 
The extended model considered ten specific points of inter- 
section and eight speed ratios. 

The second method may be classed as a Monte Carlo method, 
at least within Brown's definition: "...any procedure which 

involves the use of statistical sampling techniques to 
approximate the solution of a ... physical problem. ..." /6/ 
The Monte Carlo method utilized virtually the same computer 
program subroutines as the numerical integration approach, 
but selected the submarine's initial positions from a 
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uniform pseudorandom number (URN) generator [21], and 
selected track intersection points using the Box-Mueller 
method of normal random number generation [28] . For this 
technique, the basic model took 1,000 samples for each of 
nine speed ratios per crossing angle, and for 17 crossing 
angles. The extended model took 10,000 samples for each 
of nine speed ratios. The total number of vector 
( "maneuvering-board" ) solutions, then, exceeded 500,000; 
the usefulness of the computer in this analysis is obvious. 

1 . The Random Encounter Model 

Consider a submarine patrolling a barrier of length 
L miles at a constant speed, v. Assume that the position 
of the submarine on its patrol line at any time, x , is a 
uniformly distributed random variable, Y^. Without loss of 
generality, then, the problem can be considered to begin at 
time x when a surface vessel on a course a, projected to 
intersect the submarine patrol line at its center, is a 
distance x^/cos a from the point of intersection. Define 
the distance x^ to be equal to the range of the submarine 
detection (assuming a definite range law) , R, which in the 
initial case was considered to have been 30 miles. Refer 
to Figure 1 and the definitions summarized below: 



u = u cos a i + u sin a j 

-v t 

v = ± v J 



- u 



ship velocity ; 

submarine velocities; 

relative velocity of the 
submarine when it is on the 
up leg; 
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Figure 1. Vector Diagram of the 
Encounter (Relative Space) 
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u 



W 2 = V 2- 



t, = (y + L - Y, )/v 
1 xjr p b ' 



ti = (V b - y p )/v 



T = L/v 



R 



1 



/v 

x i + 

P 



j 




R^ + w^t 



= relative velocity of the 
submarine when it is on the 
down leg; 

= time remaining for the 
submarine to be on its 
initial leg before course 
reversal ; 

= time remaining for the sub- 
marine to be on its initial 
downleg before course 
reversal ; 

= length of time required for 
the submarine to travel 
entire length of its patrol 
line ; 

= initial position vector of 
the submarine from the ship; 

= position of the submarine 
when it makes its course 
reversal . 



2 . Algorithm to Determine Probability of Encounter 

(1) Select the initial value of Y, such that Y, = y . 

b b J p 

(2) Specify the initial direction of submarine velocity; 
call it v. 

(3) Compute the values of t, t, , R^ , and R; set w = w^. 

(4) Compute the dot product R^*w = z: 

(a) If z >.0, CPA = R^ ; go to step (5) ; 

(b) If z < 0, compute the dot product R n *w = z'; 



1) 


If z' 


•* 

o 

Vl 


CPA = R ; 

n 


2) 


If z' 


> 0, 


the CPA lies between R, and R , and 

1 n 



must be determined; 

a) set a = 1/2; 

b) compute y = (aw)/|w|; 

c) compute the dot product (R^ + y) *w = z"; 
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1. ) if z" <0, add 1/2 to a, and go to 

step 2)b) above; 

2. ) if z" >_ 0 , then CPA = R^ + y; go to 

step (5) . 

(5) (a) If | CPA | < R, a detection occurred; go to step 

(1) if the number of samples taken is greater than N, 
otherwise go to step (6). 

(b) If | CPA | < R, there was no detection on this leg: 

1) if w = w^, set w = W£ ; 

_ . -> -> -> -* 

2) if w = W£ , set w = w-^; 

3 ) go to step (4 ) . 

(6) P (Encounter) = # detections/ N. 

The general flow diagrams describing the logic for the 
models are contained in Appendix B. Note that there are 
two models, namely, the Basic Model, in which the ship track 
crosses the midpoint of the submarine patrol line at various 
angles, and the Extended Model, in which the crossing angle 
is fixed at 90 degrees, but the crossing position is varied 
along the track (and its extension) . Appendix C contains 
the actual computer instructions used to generate the data 
analyzed herein. 

D. DESCRIPTION OF THE COMPUTER SYSTEM 

The computer system used in the development of this 
thesis was the high-performance IBM System/360 Model 67, 
which is located at Ingersoll Hall, U.S. Naval Postgraduate 
School, Monterey, California. This sophisticated data 
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processing system employs a FORTRAN G compiler, along with 
the standard software and hardware associated with the 360 
system. 

Programming was accomplished using the standard FORTRAN 
IV language as proposed by the American Standards Associa- 
tion X3.4.3 FORTRAN Working Group [29]. 

Approximate execution times for the four programs in- 
cluded in this thesis were as follows: 

Basic Model: 45 minutes; 

Extended Model: 30 minutes; 

Random Sampling, Midpoint: 110 minutes; 

Random Sampling, Normally-distributed crossing: 

46 minutes. 
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II. THE RANDOM ENCOUNTER - ANALYSIS 
BY NUMERICAL INTEGRATION - 

A. RANDOM ENCOUNTER MODEL (FIXED CROSSING POINT) 

1 . Description of the Model 

The algorithm described in section l.B. was used, 
except that the values of Y^ , and the initial direction of 
the submarine were determined using numerical integration 
techniques as follows: 

The submarine's initial position for the first 
sample was established at the bottom of the patrol line, 
with the initial direction "up." The relative motion prob- 
lem was solved numerically to determine if a detection 
occurred for this configuration. For each subsequent 
sample, the submarine's initial position was incremented 
by 500 yards (1/4^^ of a nautical mile) , until each 500 
yard increment in the initial direction had been sampled, 
up to the initial position lying on the patrol line extrem- 
ity. The submarine direction was then reversed and the 
procedure continued as before except that the succeeding 
initial positions resulted from 500-yard increments in the 
new direction. The number of samples, then, was 8 x (length 
of the submarine patrol line in nautical miles). The basic 
parameters of the model from which most of the data in 
Appendices C and D were obtained were 

submarine patrol line length: 120 miles 

submarine patrol speed: 10 knots 

definite range law detection range: 30 miles 
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These parameters were used in both the basic and extended 
models for the simulation by numerical integration, and for 
the Monte Carlo simulations. 

2 . Results 

The data resulting from the simulation by numerical 
integration are contained in Appendix A; these data have 
been plotted in Figure 2 and Figure 3. Figure 2a shows how 
the detection probability varied with the ship-submarine 
speed ratio for a ship crossing the patrol center at various 
angles. The points at which detection probability was 
reduced by 90% and 95% have been connected by the two curves 
labelled "90%" and "95%," respectively.* These curves are 
shown in Figure 2b. Figure 3 is an expansion in scale of 
Figure 2b, showing greater detail. Figure 4 displays a plot 
of the Probability of Detection vs. Track Crossing Angle for 
a midpoint crossing. The behavior of P (D) as a function of 
a is of little practical interest for the case where a 
submarine commander has freedom to select a patrol line 
orientation with respect to a convoy route; a patrol 



The results plotted in Figure 2 seems to indicate that 
a barrier unit's best tactic would be to make the patrol line 
cross the ship's track at a shallow angle. This would 
indeed be the case if the ship's track were known with cer- 
tainty, as was assumed in the basic model discussed here. 

The more realistic assumption of uncertainty in the know- 
ledge of the ship's track makes the choice of a shallow 
angle undesirable from the viewpoint of the barrier unit. 

The design of an "optimal" barrier patrol line, the crossing 
angle of which depends on, say, the probability distribu- 
tion of a convoy route displacement, appears to be an 
interesting problem, but is not considered herein. 
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Figure 2a. Probability of Random Encounter vs. Ship- 
Submarine Speed Ratio (Midpoint Crossing) 



Submarine Speed: 10 knots 

Patrol length: 120 miles 

a = Crossing angle in degr 
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Probability of Detection P(D) 



Submarine Speed: 10 kts. 

Patrol Length: 120 miles 

a = Crossing angle in degrees 




Figure 3. Probability of Random Encounter vs. 
Ship-Submarine Speed Ratio (Midpoint 
Crossing) Enlarged Scale 
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perpendicular to the expected surface ship track would 
normally be selected. Figure 4 could be useful if the 
submarine were restricted for some reason to orient its 
patrol line at an angle other than 90 degrees to the convoy 
route. The effect of relaxing the assumption that the 
surface ship track crosses the patrol line at the midpoint 
of the latter is examined in the next section. The assump- 
tion of a submarine patrol line length of 120 miles (four 
times the detection range) agrees with Neuter' s "total 
attack width" which is derived in his paper. /!/ 



B. EXTENSION OF THE MODEL TO VARIABLE CROSSING POINTS 
1 . Description of the Model 

The basic model contains the expression 

y = x tan a - L/2 
P 

to insure that the ship's track will bisect the submarine's 

patrol line. It did not seem unreasonable to assume that a 

submarine normally attempts to center its patrol on the 

expected convoy route; that is, the patrol is centered on 

an expected value. A normal distribution of tracks, each 

perpendicular to the submarine patrol line, was assumed for 

the extension of the model, so that Y is a normally- 

P 

distributed random variable. In terms of the coordinate 
system used herein, 



E(Y ) = -L/2 
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The logic used in the computer program is illustrated in 
Appendix B. The results of the computer calculations are 
plotted in Figure 5. 

2 . Analysis of Results 

The flattening of the curves for the speed ratios 
p = 1 and p = 9 suggested the possibility of two simple 
limiting shapes, namely the shapes resulting when p = 0 
and p = °°. As p 0, i.e., when the ship's speed u is such 
that u<<v, it can be seen that the probability of detection 
will approach unity. More precisely, for a given detection 
radius, submarine speed, and submarine patrol line length, 
there exists some value of ship speed such that the sub- 
marine will always have an opportunity to close within 
detection range. Let 

L = length of the submarine patrol line 

y = y-coordinate of the bottom limit of the sub- 

marine patrol line 

u = ship speed 

v = submarine speed 

p = ship-submarine speed ratio 

r = detection radius (definite range law) 

P (D) = probability of detection 

Then, for all u, such that 

2r/u £ 2L/v , 

or, for all p*, such that 



r/L £ p* 
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Figure 5. Probability of Detection vs. Crossing 
Position (Perpendicular Tracks) 
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lim P (D) = 1 



v * 

p p 

if I y p I < L + r . 

For the case of p -*■ 00 , i.e., submarine speed approaching 
zero, the probability of detection is simply the probability 
that the submarine's position (which has been assumed to be 
uniformly-distributed along the patrol line) falls within 
distance r of the ship's track. Refer to Figure 6 below. 






P (D) = (r + Ay)/L 



P (D) = (r - Ay)/L 



P (D) = 2r/L 



L 






Figure 6. Probability of Detection vs. Crossing 
Position for Zero Submarine Speed 



It can be seen from Figure 6 that if the ship's track 
crosses the patrol line at a distance from the patrol line 
end-points greater than or equal to r, then the detection 
probability is simply 

2r/L 
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If the ship's track crosses the patrol line at a distance 
less than r, say Ay, from a patrol line end-point, then the 
detection probability is given by 

(r + Ay)/L 

Finally, if the ship's track does not cross the patrol line, 
but passes within detection range r, say Ay, of one of its 
end-points, the detection probability is given by 

(r - Ay)/L 

Using the convention that the coordinate system origin lies 
on the ship's position at the instant that the ship becomes 
"detectable," the above argument may be formulated as 
follows : 




The limiting case result, as presented herein, appears at 
first glance to disagree with the limiting case result ob- 
tained by Koopman ; Koopman ' s result for a zero-speed 
barrier unit is P (D) = 0.333... /8/, and the result from 
the geometric integration for p = 00 is 0.4 88. This latter 
value is a consequence of the assumption of normality; if 
Koopman ' s underlying assumption, namely a uniform 
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distribution of crossing points, were to be applied to this 
model, the limiting probability of detection would be 
(area under curve for p = °°)x (density for uniform distri- 
bution) or, 

(0.5) (60) , (0.5) (60) _ 60 _ n 

180 180 180 * •*' ' 

which is identically Koopman ' s result. The model expressed 
by equations (1) is, therefore, consistent with Koopman. 

A simple FORTRAN computer program was written (Appendix F) 
to compute the detection probability using equations (1) ; 
the results are plotted in Figure 7. Figure 7 also contains 
a plot of the computed results for comparison with the 
limiting cases. The agreement with equations (1) is evident. 
The behavior of Probability of Detection vs. Crossing 
Position in the region of ship-submarine speed ratios 
greater than three is unexpected (one would expect the 
curves to remain uninodal) . At values of p such that p>4, 
the probability of detection resulting from a patrol, the 
center of which was offset from the expected track crossing 
by about 1/6 ^ of the length of the patrol line, resulted in 
a higher detection probability than when the patrol was 
centered on the track crossing point. There data indicate 
that a submarine may effectively reduce the speed ratio by 
properly positioning itself. 
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TABLE I 



COMPARISON OF PEAKS ON PROBABILITY CURVES 
FOR SHIP-SUBMARINE SPEED RATIOS OF 5 AND 9 



SPEED 

RATIO 


P (D) FOR CENTERED 
PATROL LINE 


P (D) FOR OFF-CENTERED 
PATROL LINE (20 MILES) 


5 


.5104 


.5604 


9 


.5021 


.5292 



Table I, extracted from the data, reveals that the prob- 
ability of detection of a ship with a nine-to-one speed 
advantage over the patrolling submarine by a submarine with 
a 20-mile offset patrol line is greater than the prob- 
ability of detection of a ship with a five-to-one speed 
advantage by the same submarine with its patrol line 
centered on the crossing track. The table shows that with 
a 20-mile offset, the submarine would enjoy a higher 
probability of detection than it would for the case of a 
target with P = 5 and no offset. Figure 8 employs a con- 
tracted horizontal scale which serves to emphasize the nodal 
behavior discussed above. The conditional probability, 

P (Detection | Y = y) , 

P 

was determined by numerical analysis, and is plotted in 
Figure 7 for six values of p. The track crossing point was 
assumed to be normally-distributed, that is, N(0,400). Now 



41 




Track Crossing Position Offset (miles) 

Figure 8. Probability of Detection vs. Track Crossing 
Position (Perpendicular Tracks) 
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P (Detection) = P (R = r) , 

where r e [0,°°), and R is the random variable representing 
the distance between the submarine and the surface ship. In 
the problem at hand, r = 30 miles, and since a definite 
range law was assumed, a detection occurs when R = 30 , hence 

P (Detection) = P(D) 

is a mass function. From Baye 1 s formula, 

f r y (r,y) = f(r|Y = y)f y (y) . (2) 

' P P 

Integrating over the range of y, 

OO 00 

/ f „ ( r , y ) dy = / f(r|Y = y)f„ (y)dy . (3) 

-oo k,x p _<» p * p 

The right side of equation (3) is the expectation of the 

conditional density f(r|Y = y) , and the left side is the 

P 

marginal density, in this case the mass density, of r; hence; 

OO 

P (D) = p (r) = / f(r Y = y)f (y)dy . (4) 

K -oo P ip 

Since the conditional density in equation (4) was derived 
numerically , and no explicit mathematical expression was 
determined, the integral indicated was performed by geo- 
metric integration techniques. Appendix E shows how the 
area under the curve 

z = f (r|Y p = y) fy (y) 
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was obtained. Equation (4) may be rewritten as 



® L/2+r 

/ f (r | Y = y ) f (y ) dy = 2 / f(r|Y =y)f (y)dy, (5) 

P P 0 P 



since 



f(r|Y = y) = 0 , \j y£[- L/2 - r, L/2 + r] , 

and the integrand is symmetric. 

Figure 9 is a plot of the function to be integrated 
the values of p used are 0, 1, 1.5, 2, 3, 5, 9, and °°. The 
results of the geometrical integration are listed below in 
Table II, and are plotted in Figure 10. 



TABLE II 



GEOMETRIC INTEGRATION RESULTS 



SHIP-SUBMARINE 
SPEED RATIO 


P (D) 


0 


.9980 


1 


.9415 


1.5 


. 7895 


2 


.6816 


3 


.5781 


5 


.5233 


9 


.5008 


OO 


.4885 
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Extended Models P (D) 



Cursory examination of Figure 10 reveals that for speed 
ratios less than four, the assumption that a fixed crossing 
point will yield more pessimistic results in terms of 
higher probabilities of detection than those arising from 
the assumption of normally-distributed cross points is not 
valid. The length of the submarine patrol line was 
varied from 120 miles down to 30 miles in order that the 
effect of changing the ratio of patrol line length to detec- 
tion range might be examined. Figure 11 shows the results 
for a patrol line length of 90 miles (three times the 
detection range) . Note that the binodal behavior of the 
probability curves which was present when a patrol line 
length of 120 miles was used is absent in this case, except 
for the situation where p = 9, where the binodal behavior 
is incipient. The nodal phenomenon observed indicates a 
possible area of interest for additional research. 
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Figure 11. Probability of Detection for Various 
Speed Ratios vs. Track Crossing Position 
(Perpendicular Tracks) for L=90 miles 
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III. THE RANDOM ENCOUNTER - SIMULATION 
BY MONTE CARLO METHODS - 

A. MIDPOINT CROSSING MODEL 

The encounter model developed by deterministic methods 
was used to build the models and computer programs for the 
stochastic simulation. All assumptions pertaining to the 
deterministic model were included in the Monte Carlo ap- 
proach. 

The submarine's initial position on his patrol line was 
determined for each sample taken by selecting from a uniform 
(0,1) pseudorandom number generator. This procedure agreed 
with the assumption that the submarine was equally likely to 
be anywhere on his patrol line at any time t. 

The submarine's initial up or down heading on its patrol 
line was determined by selecting a pseudorandom number from 
a uniform (0,1) variate generator. If this number was less 
than 0.5 the submarine was assigned an initial up heading; 
if the number was greater than 0.5 the submarine was 
assigned an initial down heading. If the random variate was 
equal to 0.5 another uniform random variable was selected 
and the above procedure was repeated. 

The midpoint crossing of the patrol line by the merchant 
was arranged by positioning the bottom of the patrol line in 
such a manner as to insure that this phenomenon occurred. 
Section (11) gives the analytical explanation dealing with 
the development of the deterministic approach. 
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The results obtained using the Monte Carlo technique are 
presented in Appendix A. Graphs depicting the results can 
be found in Figure 12. It should be noted that the results 
obtained utilizing the stochastic simulation agree with 
those obtained from the deterministic approach. 

B. EXTENSION TO NORMALLY-DISTRIBUTED CROSSING POINTS 

The model and computer program for normally-distributed 
crossing points of the convoy route about the center of the 
submarine patrol line is an extension of the midpoint 
crossing model. 

The submarine's initial position and heading were 
determined as for the midpoint crossing model. The point 
of intersection of the submarine patrol line and the convoy 
route was determined by generating a normal N( 0,400) 
pseudorandom variate using the Box-Mueller Method [28]. The 
mean of the normal variate corresponded to the center of the 
patrol line. The problem was then generated in a manner 
similar to that of the midpoint crossing model. 

A relatively small sample size of 10 normal variates, 
one for each 1,000 uniformly-distributed initial starting 
positions for the submarine, was generated due to the time 
limitations on the computer system. This sample size took 
approximately 46 minutes to execute. The computer program 
which appears in Appendix D has been written using a larger 
sample size in an effort to aid in any future duplications 
of the simulation. 
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Speed Ratio p=u/v 

Figure 12. Detection Probability vs. Ship-Sub Speed Ratio 
with Parameter of Route Angle (Monte Carlo technique) 



The results obtained using this simulation technique 
are tabulated in Appendix A. A plot of the results appears 
in Figure 13 and Figure 14. The results and analysis of 
this situation are similar to those obtained from the method 
of analysis by numerical integration in the extension to 
variable crossing points as is evident from Figure 14. 
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Sub Speed Ratio as Parameter (Monte Carlo Technique) 
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Figure 14. Detection Probability vs. Ship-Sub Speed Ratio for 
Route Angle a=0° for Midpoint Crossing and Normalized 
Convoy Route Approaches 



IV. AN EXAMPLE OF A NORTH ATLANTIC CROSSING 



A. ASSUMPTIONS 

It was assumed that the major western North Atlantic 
seaports are New York, Norfolk, Windward Passage, and 
Puerto Rico (Panama Canal traffic) ; the major eastern 
North Atlantic seaports were assumed to be the Faero Islands 
(representing the straits between Iceland and Scotland) , 
the English Channel, and Gibraltar. This assumption results 
in twelve possible Atlantic-crossing convoy routes. These 
possible convoy routes were used as a basis to determine 
the following assumptions: 

1) allocation of enemy submarines 

2) lengths of submarine barrier patrol lines 

3) angle of orientation of the submarine patrol lines 
to the convoy tracks . 

A consequence resulting from the assumption of specific 
convoy routes is a deviation from the assumptions made by 
Koopman [18], Neufer [5], Dobbie [24], and others which was 
that of a uniform distribution of submarines throughout an 
area. Figure 15 is a sketch which illustrates the rationale 
for the assumptions made herein. 

Given that the only possible convoy routes are as shown 
in the figure, an enemy might be expected to allocate his 
anti-shipping submarine patrols as follows: 

1. One submarine to patrol each port (since ASW 
activity would be expected to be heaviest near these ports. 
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greater numbers of patrols at these locations was not 
assumed) . 

2. A submarine patrol at each expected juncture of 
convoy routes (this effectively places a barrier on each 
of the convoy routes which merge at the junction) . 

3. No segment of convoy track greater than 500 miles 
to be unopposed by a barrier. 

The preceding assumptions require that an enemy have 
the capability to maintain 70 submarines on station in 
the North Atlantic. Note that this distribution results 
in 105 barriers against the possible convoy routes assumed, 
even though only 70 submarines are employed. 

Another assumption made was that the submarine barriers 
would be oriented according to the following rules: 

1. Patrolling of straits to be oriented so as to 
minimize the patrol line length. 

2. Patrolling of seaport exits (entrances) to be per- 
pendicular to, and centered on, a bisector of the widest 
angle formed by two departure courses. 

3. Patrolling of route junctions to be along the 
bisector of the largest angle formed by the intersection 
of the two convoy routes, and to be centered on the 
expected junction. 

4. All other patrols to be centered on, and oriented 
at right angles to the expected route. 

Additionally, all submarine patrolling speeds were assumed 
to be 10 knots, and all patrol lines were assumed to be 
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120 miles in length (unless shortened by geographical 
limitations ) . 

B. APPLICATION OF THE MODELS 

1 . Orientation of Barriers on the Norfolk-English 
Channel Route 

The route between Norfolk and the English Channel 
is approximately 3,500 miles long, and intersects the routes 
from New York to Gibraltar, Windward Passage to the Faero 
Islands, and Puerto Rico to the Faero Islands. The rules 
assumed required that the barrier patrols on this convoy 
route be given the following parameters: 



TABLE III 

PARAMETERS OF THE BARRIERS FOR THE EXAMPLE 



BARRIER 

# 


BARRIER LOCATION 


a° 


L 


1 


Norfolk 


5 


120 miles 


2 


N.Y. -Gibraltar junction 


12 


120 


3,4 


Along the route 


0 


120 


5 


W.P. -Faero junction 


10 


120 


6 


P.R. -Faero junction 


15 


120 


7 


English Channel 
(Land's End-Ushant Is.) 


0 


90 



2 . Solution to the Example Problem 
In general, 

P (Detection) = 1 - P (No Detection) , 
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and in this case, 

P (No Detection) = P (No Detection by Barrier #1) 

xP (No Detection by Barrier #2) 
x ... x 

xP (No Detection by Barrier n) , 
assuming stochastic independence for simplicity. In this 
example, n = 7. Since the probability of detection is a 
function of the angle of track a, the patrol line length L, 
the detection range r, and the ship-submarine speed ratio 
p, the probability of detection for barrier i may be 
represented as 

P i (D) = P(D;a i ,L i ,r i ,p i ) , 

and 

7 

P (D) = n [1 - P(D;a i ,L i ,r i ,p i ) ] . 

i=l 

According to the assumptions made. 



r 



1 



r~ = 



r 



P ^ P2 ... P ^ p , 



and, 

7 

P (Detection) = II [1 - P (D; a. ,L. ,r , p) ] 

i=l 1 1 



( 6 ) 



The probabilities of detection for each of the barriers 
listed in Table 3 are tabulated in Table IV. 
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TABLE IV 
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Using the results from Table IV in equation (6) , 



1 - P(D|p=p*) = [1 - P(D|p=p*,L = 120 , a = 0 , r = 30)] 6 

x[l - P (D | p=p* ,L = 90,a=0,r = 30)] 

The results are tabulated in Table 5. 



TABLE V 

P (DETECTION) VERSUS SPEED FOR 
THE ROUTE OF THE EXAMPLE 



p * 


P (D | p=p*) 


0 


1.0000 


1 


1.0000 


1.5 


1.0000 


2 


1.0000 


2.5 


1.0000 


3 


.9996 


4 


.9986 


5 


.9978 


9 


.9957 


oo 


.9910 



Table V indicates, ceterus paribus, that speed advantage 
alone for the surface ship will not lessen the probability 
of detection significantly by patrolling submarines for a 
typical North Atlantic crossing. 
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V. SUMMARY AND CONCLUSIONS 



A. MONTE CARLO VERSUS NUMERICAL INTEGRATION METHODS 

Figure 14 shows that close agreement exists between the 

Monte Carlo and numerical integration techniques; the numer- 
ical integration method is more economical in terms of 
computer usage time. 

B. THE INDIVIDUAL RANDOM ENCOUNTER 

1 . Midpoint Crossing at Various Angles 

Detection probability decreased from unity at P i 1 
for all angles of intersection between the ship's track and 
the submarine barrier patrol line, to 0.5 at p > 9 (for 
intersection at right angles) . Ninety-five per cent of the 
decrease occurred at speed ratios of eight and below for 
all angles, and for speed ratios of 4.5 and below for inter- 
sections at right angles. Stated another way, for assumed 
submarine speed of 10 knots, and using detection probability 
as the sole criterion for selecting ship speed capabilities, 
merchant ship speeds greater than 40 knots do not seem 
warranted. 

2 . Normally-Distributed Crossing Points 

a. The assumption of normally-distributed crossing 
points resulted in lower detection probabilities for ship- 
submarine speed ratios less than four, than did the assump- 
tion of a fixed crossing point at the center of the patrol 
line. For speed ratios greater than four, both assumptions 
yielded essentially identical results. 
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b. If ship-submarine speed ratios equal to or 
greater than four are to be considered, then the binodal 
behavior of the probabilities of detection, illustrated in 
Figure 8, should be examined further. It appears that 
this behavior may also be a function of the ratio of sub- 
marine barrier patrol length to assumed detection range. 

C. NORTH ATLANTIC TRANSIT EXAMPLE 

When six or more submarine barriers are crossed by the 
merchant ship's track, the reduction in P (Detection) 
resulting from increasing the merchant ship's speed is 
negligible. It is evident that, from the viewpoint of 
decreasing random detection probability, a large speed 
advantage is not the answer. If a large speed advantage 
by a merchant ship over a patrolling submarine increases 
the survivability of that merchant ship, then the increase 
in survivability is due, not to decreased detections, but 
to other factors, e.g., reduced submarine weapon system 
effectiveness or reduced exposure time. 
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APPENDIX A 



PPCBABIL ITY CF R ANDCM ENCCUNTER 
BASIC MODEL 

PlFINITc. RANGE LAW 

CFTfCTICN RANGE: 30. C MILES 

LENGTH OF SLB PATRCL LINE: 120.0 MILES 

X-C' I STANCE TO PATRCL LINE: 30.0 MILES 

SUBMARINE SPEED: 1G.0 KNCTS 

NO. OF SAMPLES TAKEN FCR EACH SCENARIO: 96C 



A VO . NO. 


SHI P 


SPEED 


RELATIVE 


SPEECS 


RCUTE 


PROBABILITY 


LEGS 


spier 


RATIO 


UPLEG 


DCWNLEG 


ANGLE 


DETECTION 


I • 4 C 


1 c.c 


1.0 


14.1 


14.1 


0 


1.0000 


1.36 


1 5 • C 


1.5 


18. C 


18.0 


0 


0.9354 


1.3? 


2 C . 0 


2.0 


22.4 


22.4 


0 


C. 8000 


1.2P 


2 5,0 


2.5 


26.9 


26. 9 


0 


0.6937 


1.25 


30. ? 


3.0 


31.6 


31.6 


0 


0.6250 


1.22 


3 5.0 


3.5 


36.4 


36.4 


0 


0.5750 


1.21 


40.0 


4.0 


41.2 


41.2 


0 


0.5437 


1.1 6 


5C.C 


5.0 


51. C 


51.0 


0 


0.5104 


1.17 


6C.0 


6.0 


6 C . 8 


60.3 


0 


0.5062 


1.15 


7C.0 


7.0 


70.7 


70.7 


0 


0.5062 


1.15 


ec.o 


e.o 


30.6 


80.6 


0 


0.5021 


1 • 1 A 


9 0.0 


9.0 


9C.6 


90.4 


0 


0.5021 


1.39 


1C.0 


1.0 


10.7 


16.9 


25 


1.0000 


1.26 


15.0 


1.5 


14.1 


21.3 


25 


0.9427 


1.32 


2C.0 


2.0 


18.2 


25.9 


25 


C. 8698 


1.29 


25.0 


2.5 


22.7 


30.6 


25 


0.7781 


1.27 


3 C . 0 


3.0 


27.3 


35.4 


25 


0.7240 


1.2A 


35.0 


3. 5 


32.1 


40.3 


25 


0.6719 


1.22 


4 C . 0 


A.O 


36.9 


45.1 


25 


0.6406 


1.19 


5C.0 


5.0 


46.7 


55.0 


25 


0.5865 


1.17 


6 C.C 


6.0 


56.5 


64.9 


25 


0.5687 


1.16 


7C.C 


7.0 


66.4 


74.8 


25 


0.5562 


1.15 


80.0 


8.0 


76.3 


84.7 


25 


0.5542 


1. 14 


9C.0 


9.0 


86.3 


94.7 


25 


0.5542 
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AVG.NO. 


SHIP 


SPEED 


RELATIVE 


SPEEDS 


ROUTE 


PROBABILITY 


LEGS 


SPEED 


RATIO 


UPLEG 


DOWNLEG 


ANGLE 


DETECTION 


1.39 


10.0 


1.0 


10. C 


17.3 


30 


1.0000 


1.36 


15.0 


1.5 


13.2 


21.8 


30 


0.9469 


1.32 


20.0 


2.0 


17.3 


26.5 


30 


0.8833 


1.28 


25.0 


2.5 


21.8 


31.2 


30 


0.8187 


1.26 


30.0 


3.0 


26.5 


36.1 


30 


0.7615 


1.25 


35.0 


3.5 


31.2 


40.9 


30 


0.7229 


1.23 


40.0 


4.0 


36.1 


45.8 


30 


0.6865 


1.20 


50.0 


5.0 


45.8 


55.7 


30 


0.6344 


1.18 


60.0 


6.0 


55.7 


65.6 


30 


0.6052 


1.16 


7C.0 


7.0 


65.6 


75.5 


30 


0.5906 


1.15 


80.0 


8.0 


75.5 


85.4 


30 


0.5802 


1.15 


90.0 


9.0 


85.4 


95.4 


30 


0.5792 


1.39 


10.0 


1.0 


9.2 


17.7 


35 


1.0000 


1.35 


15.0 


1.5 


12.4 


22.3 


35 


0.9552 


1.31 


20.0 


2.0 


16.4 


27.0 


35 


0.8969 


1.28 


25.0 


2.5 


20.9 


31.8 


35 


0.8583 


1.26 


30.0 


3.0 


25.6 


36.7 


35 


0.8083 


1.24 


35.0 


3.5 


30.4 


41.6 


35 


0.7677 


1.23 


40.0 


4.0 


35.2 


46.5 


35 


0.7396 


1.20 


50.0 


5.0 


45.0 


56.3 


35 


0.6948 


1.18 


60.0 


6.0 


54.9 


66.2 


35 


0.6583 


1.17 


70.0 


7.0 


64.8 


76.2 


35 


0.6375 


1.16 


80.0 


8.0 


74.7 


86.1 


35 


0.6250 


1.15 


90.0 


9.0 


84.7 


96.1 


35 


0.6177 


1.39 


10.0 


1.0 


8.5 


18.1 


40 


1.0000 


1.34 


15.0 


1.5 


11.5 


22.8 


40 


0.9677 


1.31 


20.0 


2.0 


15.6 


27.5 


40 


0.9177 


1.28 


25.0 


2.5 


20.1 


32.3 


40 


0.8812 


1.25 


30.0 


3.0 


24.8 


37.2 


40 


0.8594 


1.23 


35.0 


3. 5 


29.6 


42.1 


40 


0.8323 


1.22 


40.0 


4.0 


34.4 


47.1 


40 


0.8021 


1.20 


50.0 


5.0 


44.2 


56.9 


40 


0.7573 


1.19 


60.0 


6.0 


54.1 


66.9 


40 


0.7323 


1.17 


70.0 


7.0 


64.0 


76.8 


40 


0.7052 


1.16 


80.0 


8.0 


74.0 


86.8 


40 


0.6906 


1.15 


90.0 


9.0 


83.9 


96.7 


40 


0.6750 
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4VG.N0. 


SHIP 


SPEED 


RELATIVE 


SPEEDS 


ROUTE 


PROBABILITY 


LEGS 


SPEED 


RATIO 


UPLEG 


DOWNLEG 


ANGLE 


DETECTION 


1.39 


10.0 


1.0 


7.7 


18.5 


45 


1.0000 


1.34 


15.0 


1.5 


10.6 


23.2 


45 


0.9865 


1.3C 


20.0 


2.0 


14.7 


28.0 


45 


0.9365 


1.27 


25.0 


2.5 


19.3 


32.9 


45 


0.9115 


1.25 


30. 0 


3.0 


24.0 


37.7 


45 


0.8927 


1.23 


35. C 


3. 5 


28.8 


42.7 


45 


0.8802 


1.21 


4C.C 


4.0 


33.7 


47.6 


45 


0.87CR 


1.19 


50.0 


5.0 


43.5 


57.5 


45 


0.8344 


1. 18 


6C.0 


6.0 


53.4 


67.4 


45 


0.8073 


1.17 


70.0 


7.0 


63.3 


77.4 


45 


0.7875 


1.16 


80.0 


8.0 


73.3 


87.4 


45 


0.7740 


1.16 


90.0 


9.0 


83.2 


97.3 


45 


0.7625 


1.39 


10.0 


1.0 


6.8 


18.8 


50 


1.0000 


1.33 


15.0 


1.5 


9.8 


23.6 


50 


l.OOCO 


1.28 


2C.0 


2.0 


13.9 


28.4 


50 


0.9656 


1.25 


25.0 


2.5 


18.5 


33.3 


50 


0.9427 


1.23 


3 C. 0 


3.0 


23.2 


38.2 


50 


0.9302 


1.21 


35.0 


3.5 


28.1 


43. 1 


50 


0.9240 


1.20 


4C.0 


4.0 


33.0 


48. 1 


50 


0.9177 


1. 19 


50. C 


5.0 


42.8 


58.0 


50 


0.9083 


1. 17 


60.0 


6.0 


52.7 


68.0 


50 


0.9021 


1. 16 


70.0 


7.0 


62.7 


77.9 


50 


0.8854 


1.15 


80. 0 


e.o 


72.6 


87.9 


50 


0.8698 


1. 14 


90. C 


9.0 


82.6 


97.9 


50 


0.8583 


1.39 


1G.C 


1.0 


6.0 


19.1 


55 


1.0000 


1.31 


15.0 


1.5 


8.9 


23.9 


55 


1.0000 


1.26 


20.0 


2.0 


13.1 


28. 9 


55 


1.0000 


1.22 


25.0 


2.5 


17.8 


33.7 


55 


0.9844 


1.20 


30.0 


3.0 


22.6 


38.6 


55 


0.9719 


1.16 


35.0 


3.5 


27.4 


43.6 


55 


0.9656 


1.17 


40. C 


4.0 


32.3 


48.5 


55 


0.9615 


1. 15 


5C.0 


5.0 


42.2 


58.5 


55 


0.9552 


1.13 


60.0 


6.0 


52.1 


68.4 


55 


0.9510 


l. 12 


70.0 


7.0 


62. 1 


78.4 


55 


0.9490 


1.12 


80.0 


8.0 


72.0 


88.4 


55 


0.9469 


1.11 


90.0 


c .0 


82.0 


98.4 


55 


0.9448 
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•• J5 



* O-J, -I 



AVG.NC. 


SHI P 


SPEED 


LEGS 


SPEED 


RATIO 


1.39 


1C.0 


1.0 


1.29 


15.0 


1.5 


1.23 


20.0 


2.0 


1.19 


25.0 


2.5 


1.16 


30.0 


3.0 


1.1^ 


35. C 


3. 5 


1.12 


40. 0 


4.0 


1.10 


50.0 


5.0 


1 • C 8 


60.0 


6. 0 


1 • C 7 


70.0 


7.0 


1.C6 


80. 0 


8.0 


1.C5 


90.0 


9.0 



RELATIVE SPEEDS 
UPLEG DOWNLEG 

5.2 19.3 

8.1 24.2 

12.4 29. 1 

17.1 34.0 

21.9 39.0 

26.8 43.9 

31.7 48.9 

41.6 58.9 

51.6 68.8 

61.5 78.8 

71.5 88.8 

81.5 98.8 



ROUTE 


PROBABILITY 


angle: 


DETECTION 


60 


1.0000 


60 


> 1.0000 


60 


1.0000 


60 


1.0000 


60 


1.0000 


60 


1.0000 


60 


1.0000 


60 


1.0000 


60 ' 


1.0000 


60 


1.0000 


60 


1.0000 


60 


1.0000 
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PROBABILITY CF RANDOM ENCOUNTER 
EXTENDED MODEL 



DEFINITE RANGE LAW 

DETECTION RANGE: 30.0 MILES 

LENGTH OF SUB PATRCL LINE: 120.0 MILES 

X-DI STANCE TO PATRCL LINE: 30.0 MILES 

CROSSING ANGLE: C.O CEGREES (PERPENDICULAR) 



SUBMARINE SPEED: 10.0 KNOTS 



NO. 


OF SAMPLES 


TAKEN FOR 


EACH SCENARIO: 


960 


AVG.NO. 


SHIP 


SPEED 


RELATIVE 


SPEEDS 


ROUTE 


PROBABILITY OF 


LEGS 


SPEED 


RATIO 


UPLEG 


DOWNLEG 


POSIT 


DETECTION 


1.40 


10.0 


1.0 


14.1 


14.1 


60.0 


l.OOOC 


1.36 


15. C 


1.5 


18.0 


18.0 


60.0 


0.9354 


1.23 


20. C 


2.0 


22.4 


22.4 


60.0 


0.8000 


1.28 


25.0 


2.5 


26.9 


26.9 


60.0 


0.6937 


1.25 


3 C. C 


3.0 


31.6 


31.6 


60.0 


0.6250 


1.21 


40. C 


*.0 


41.2 


41.2 


60.0 


0.5437 


1.18 


50.0 


5.0 


51.0 


51.0 


60.0 


0.5104 


2.32 


9C.0 


9.0 


90.6 


90.6 


60.0 


0.5021 


1.4C 


1C.G 


1.0 


14.1 


14.1 


65.0 


1.0000 


1.36 


15.0 


1.5 


18.0 


18.0 


65.0 


0. 9260 


1.32 


20. 0 


2.0 


22.4 


22.4 


65.0 


0.7844 


1.28 


25.0 


2.5 


26.9 


26.9 


65.0 


0.6937 


1.25 


30. C 


3.0 


31.6 


31.6 


65.0 


0.6250 


1.21 


^C.C 


4.0 


41.2 


41.2 


65.0 


0.5500 


1.18 


5C.0 


5.0 


51.0 


51.0 


65. D 


0.5260 


2.32 


90.0 


9.0 


90.6 


90.6 


65.0 


0.5021 


1.40 


1C.0 


1.0 


14.1 


14.1 


70.0 


l.OOOC 


1.35 


15.0 


1.5 


18.0 


18.0 


7C.0 


0.8687 


1.30 


2C.0 


2.0 


22.4 


22.4 


70.0 


0.7635 


1.26 


25.0 


2.5 


26.9 


26.9 


70.0 


0.6750 


1.24 


30.0 


3.0 


31.6 


31.6 


70.0 


0.6177 


1.21 


40.0 


4.0 


41.2 


41.2 


70.0 


0.5708 


1.18 


50.0 


5.0 


51.0 


51.0 


70.0 


0.5469 


2.32 


5C.C 


9.0 


90.6 


90.6 


70.0 


0.5021 


1.40 


10.0 


1.0 


14.1 


14.1 


75.0 


0.9781 


1.33 


15.0 


1.5 


18.0 


18.0 


75.0 


0.8062 


1.28 


20. 0 


2.0 


22.4 


22.4 


75.0 


0.7125 


1.24 


25. C 


2.5 


26.9 


26.9 


75.0 


0.6542 


1.22 


3C.0 


3.0 


31.6 


31.6 


75.0 


0.6104 


1.19 


40.0 


4.0 


41.2 


41.2 


75.0 


0.5771 


1.18 


50.0 


5.C 


51.0 


51.0 


75.0 


0.5604 


2.31 


9C.0 


9.0 


90.6 


90.6 


75.0 


0.5229 
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avg . no . 


SHIP 


SPEED 


RELATIVE 


SPEEDS 


RCUTE 


PROBABILITY i 


LEGS 


SPEED 


RATIO 


UPLEG 


DCWNLEG 


POSIT 


DETECTION 


1.40 


10.0 


1.0 


14.1 


14.1 


80.0 


0.9365 


1.31 


15.0 


1.5 


18.0 


18.0 


80.0 


0.7437 


1.26 


20. 0 


2.0 


22.4 


22.4 


80.0 


0.6500 


1.22 


25.0 


2.5 


26.9 


26.9 


80.0 


0.6031 


1.20 


30.0 


3.0 


31.6 


31.6 


80.0 


0.5969 


1.17 


40.0 


4.0 


41.2 


41.2 


80.0 


0.5771 


1.16 


50 . C 


5.0 


51.0 


51.0 


80.0 


0.5604 


2.28 


90.0 


9.0 


90.6 


90.6 


80.0 


0.5292 


1.39 


10.0 


1.0 


14.1 


14.1 


85.0 


0.8844 


1.29 


15.0 


1.5 


18.0 


18.0 


85.0 


0.6812 


1.23 


20.0 


2.0 


22.4 


22.4 


85.0 


0.5875 


1.20 


25.0 


2.5 


26.9 


26.9 


85.0 


0.5615 


1.18 


30.0 


3.0 


31.6 


31.6 


85.0 


0.5552 


1.15 


40.0 


4.0 


41.2 


41.2 


85.0 


0.5490 


1.13 


50.0 


5.0 


51.0 


51.0 


85.0 


0.5469 


2.24 


90.0 


9.0 


90.6 


90.6 


85.0 


0.5292 


1.36 


10.0 


1.0 


14.1 


14.1 


90.0 


0.8219 


1.27 


15.0 


1.5 


18.0 


18.0 


90.0 


0.6187 


1.21 


20.0 


2.0 


22.4 


22.4 


90.0 


0.5302 


1.18 


25.0 


2.5 


26.9 


26.9 


. 90.0 


0.5198 


1.16 


30.0 


3.0 


31.6 


31.6 


90.0 


0.5135 


1.13 


40.0 


4.0 


41.2 


41.2 


90.0 


0.5073 


1.11 


50.0 


5.0 


51.0 


51.0 


90.0 


0.5052 


2.19 


90.0 


9.0 


90.6 


90.6 


90.0 


0.5010 


1.30 


10.0 


1.0 


14.1 


14.1 


110.0 


0.5719 


1.21 


15.0 


1.5 


18.0 


18.0 


110.0 


0.3844 


1.17 


20.0 


2.0 


22.4 


22.4 


110.0 


0.3635 


1.14 


25.0 


2.5 


26.9 


26.9 


110.0 


0.3531 


1.13 


30.0 


3.0 


31.6 


31.6 


110.0 


0.3469 


1.11 


40.0 


4.0 


41.2 


41.2 


110.0 


0.3406 


1.09 


50.0 


5.0 


51.0 


51.0 


110. 0 


0.3385 


2.16 


90.0 


9.0 


90.6 


90.6 


110.0 


0.3344 


1.33 


10.0 


1.0 


14.1 


14 . 1 


130.0 


0.3219 


1.23 


15.0 


1.5 


18.0 


18.0 


130.0 


0.2177 


1.18 


20.0 


2.0 


22.4 


22.4 


130.0 


0.1969 


1.14 


25.0 


2.5 


26.9 


26.9 


130.0 


0.1865 


1.13 


30.0 


3.0 


31.6 


31.6 


130.0 


0.1802 


1.11 


40.0 


4.0 


41.2 


41.2 


130.0 


0.1740 


1.09 


50.0 


5.0 


51.0 


51.0 


130.0 


0.1719 


2.16 


90.0 


9.0 


90.6 


90.6 


130.0 


0.1677 


1.58 


1C.0 


1.0 


14.1 


14.1 


150.0 


0.0 


1.42 


15.0 


1.5 


18.0 


18.0 


150.0 


0.0 


1.33 


20.0 


2.0 


22.4 


22.4 


150.0 


0.0010 


1.28 


25.0 


2.5 


26.9 


26.9 


150.0 


0.0010 


1.25 


30.0 


3.0 


31.6 


31.6 


150.0 


0.0010 


1.21 


40.0 


4.0 


41.2 


41.2 


150.0 


0.0010 


1.18 


50.0 


5.0 


51.0 


51.0 


150.0 


0.0010 


2.32 


9C.0 


9.0 


90.6 


90.6 


150.0 


0.0 
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MONTE CARLO TECHNIQUE 
PROBABILITY OF RANDOM ENCOUNTER 
MID-POINT CROSSING 

DEFINITE RANGE LAW 

DETECTION RANGE: 30.0 MILES 

LENGTH OF SUB PATROL LINE: 120.0 MILES 

X-DISTANCE TO PATROL LINE: 30.0 MILES 

SUBMARINE SPEED: 10.0 KNOTS 

NO. OF SAMPLES TAKEN FOR EACH SCENARIO: 1000 



G. NO. 


SHIP 


LEGS 


SPEED 


1.41 


10.0 


1.34 


20.0 


1.25 


30.0 


1.21 


40.0 


1.18 


50.0 


1.17 


60.0 


1.15 


70.0 


1.16 


80.0 


1.13 


90.0 


1.42 


10.0 


1.32 


20.0 


1.25 


30.0 


1.19 


40.0 


1.18 


50.0 


1.19 


60.0 


1.15 


70.0 


1.14 


80.0 


1.14 


90.0 


1.38 


10.0 


1.34 


20.0 


1.24 


30.0 


1.23 


40.0 


1.17 


50.0 


1.15 


60.0 


1.15 


70.0 


1.14 


80.0 


1.14 


90.0 



SPEED RELAT 
RATIO UPLEG 

1.0 14.1 

2.0 22.4 

3.0 31.6 

4.0 41.2 

5.0 51.0 

6.0 60.8 

7.0 70.7 

8.0 80.6 

9.0 90.6 

1.0 13.5 

2.0 21.6 

3.0 30.8 

4.0 40.4 

5.0 50.1 

6.0 60.0 

7.0 69.8 

8.0 79.8 

9.0 89.7 

1.0 12.9 

2.0 20.7 

3.0 29.9 

4.0 39.5 

5.0 49.3 

6.0 59.1 

7.0 69.0 

8.0 78.9 

9.0 88.8 



SPEED ALFA 



DOWNLEG 




14.1 


0.0 


22.4 


0.0 


31.6 


0.0 


41.2 


0.0 


51.0 


0.0 


60.8 


0.0 


70.7 


0.0 


80.6 


0.0 


90.6 


0.0 


14.7 


0.0873 


23.1 


0.0873 


32.4 


0.0873 


42.1 


0.0873 


51.8 


0.0873 


61.7 


0.0873 


71.6 


0.0873 


81.5 


0.0873 


91.4 


0.0873 


15.3 


0.1745 


23.9 


0.1745 


33.2 


0.1745 


42.9 


0.1745 


52.7 


0.1745 


62.5 


0.1745 


72.4 


0.1745 


82.3 


0.1745 


92.3 


0.1745 



PROB. 

DETECTION 

1.0000 

0.8200 

0.6270 

0.5470 

0.5290 

0.5060 

0.4810 

0.4810 

0.5070 

1.0000 

0.7650 

0.6420 

0.5460 

0.5210 

0.5020 

0.4980 

0.5180 

0.4930 

1.0000 

0.8140 

0.6380 

0.5640 

0.5190 

0.5130 

0.5210 

0.5070 

0.5000 
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f G . NO . 


SHIP 


SPEED 


RELATIVE SPEED 


ALFA 


PROB . 


LEGS 


SPEED 


RATIO 


UPLEG 


DOWNLEG 

i 




DETECTION 


1.40 


10.0 


1.0 


12.2 


15.9 


0.2618 


1.0000 


1.32 


20.0 


2.0 


19.9 


24.6 


0.2618 


0.8120 


1.26 


30.0 


3.0 


29.1 


34.0 


0.2618 


0.6500 


1.22 


40.0 


4.0 


38.6 


43.7 


0.2618 


0.5530 


1.20 


50.0 


5.0 


48.4 


53.5 


0.2618 


0.5260 


1.17 


60.0 


6.0 


58.2 


63.3 


0.2618 


0.5270 


1.15 


70.0 


7.0 


68.1 


73.2 


0.2618 


0.5190 


1.14 


80.0 


8.0 


78.0 


83.2 


0.2618 


0.5150 


1.15 


90.0 


9.0 


87.9 


93.1 


0.2618 


0.5190 


1.39 


10.0 


1.0 


11.5 


16.4 


0.3491 


1.0000 


1.32 


20.0 


2.0 


19.1 


25.2 


0.3491 


0.8600 


1.24 


30.0 


3.0 


28.2 


34.7 


0.3491 


0.6920 


1.20 


40.0 


4.0 


37.8 


44.4 


0.3491 


0.6070 


1.19 


50.0 


5.0 


47.5 


54.2 


0.3491 


0.5470 


1.18 


60.0 


6.0 


57.4 


64.1 


0.3491 


0.5480 


1.16 


70.0 


7.0 


67.2 


74.0 


0.3491 


0.5290 


1.14 


80.0 


8.0 


77.2 


83.9 


0.3491 


0.5340 


1.16 


90.0 


9.0 


87.1 


93.9 


0.3491 


0.5070 


1.42 


10.0 


1.0 


10.7 


16.9 


0.4363 


1.0000 


1.30 


20.0 


2.0 


18.2 


25.9 


0.4363 


0.8570 


1.26 


30.0 


3.0 


27.3 


35.4 


0.4363 


0.7380 


1.22 


40.0 


4.0 


36.9 


45.1 


0.4363 


0.6430 


1.18 


50.0 


5.0 


46.7 


55.0 


0.4363 


0.5910 


1.18 


60.0 


6.0 


56.5 


64.9 


0.4363 


0.5530 


1.15 


70.0 


7.0 


66.4 


74.8 


0.4363 


0.5860 


1.16 


80.0 


8.0 


76.3 


84.7 


0.4363 


0.5360 


1.13 


90.0 


9.0 


86.3 


94.7 


0.4363 


0.5490 


1.41 


10.0 


1.0 


10.0 


17.3 


0.5236 


1.0000 


1.30 


20.0 


2.0 


17.3 


26.5 


0.5236 


0.8890 


1.26 


30.0 


3.0 


26.5 


36.1 


0.5236 


0.7270 


1.24 


40.0 


4.0 


36.1 


45.8 


0.5236 


0.6690 


1.21 


50.0 


5.0 


45.8 


55.7 


0.5236 


0.6100 


1.18 


60.0 


6.0 


55.7 


65.6 


0.5236 


0.5830 


1.18 


70.0 


7.0 


65.6 


75.5 


0.5236 


0.5910 


1.15 


80.0 


8.0 


75.5 


85.4 


0.5236 


0.5780 


1.13 


90.0 


9.0 


85.4 


95.4 


0.5236 


0.5540 


1.37 


10.0 


1.0 


9.2 


17.7 


0.6109 


1.0000 


1.31 


20.0 


2.0 


16.4 


27.0 


0.6109 


0.9010 


1.26 


30.0 


3.0 


25.6 


36.7 


0.6109 


0.8140 


1.25 


40.0 


4.0 


35.2 


46.5 


0.6109 


0.7480 


1.19 


50.0 


5.0 


45.0 


56.3 


0.6109 


0.7090 


1.19 


60.0 


6.0 


54.9 


66.2 


0.6109 


0.6540 


1.17 


70.0 


7.0 


64.8 


76.2 


0.6109 


0.6190 


1.16 


80.0 


8.0 


74.7 


86.1 


0.6109 


0.6290 


1.14 


90.0 


9.0 


84.7 


96.1 


0.6109 


0.6310 
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' G . NO . 


SHIP 


SPEED 


LEGS 


SPEED 


RATIO 


1.40 


10.0 


1.0 


1.31 


20.0 


2.0 


1.24 


30.0 


3.0 


1.23 


40.0 


4.0 


1.22 


50.0 


5.0 


1.19 


60.0 


6.0 


1.16 


70.0 


7.0 


1.16 


80.0 


8.0 


1.14 


90.0 


9.0 


1.42 


10.0 


1.0 


1.30 


20.0 


2.0 


1.25 


30.0 


3.0 


1.20 


40.0 


4.0 


1.20 


50.0 


5.0 


1.20 


60.0 


6.0 


1.16 


70.0 


7.0 


1.16 


80.0 


8.0 


1.17 


90.0 


9.0 


1.38 


10.0 


1.0 


1.26 


20.0 


2.0 


1.24 


30.0 


3.0 


1.22 


40.0 


4.0 


1.17 


50.0 


5.0 


1.17 


60.0 


6.0 


1.16 


70.0 


7.0 


1.14 


80.0 


8.0 


1.14 


90.0 


9.0 


1.39 


10.0 


1.0 


1.25 


20.0 


2.0 


1.18 


30.0 


3.0 


1.16 


40.0 


4.0 


1.16 


50.0 


5.0 


1.12 


60.0 


6.0 


1.12 


70.0 


7.0 


1.11 


80.0 


8.0 


1.11 


90.0 


9.0 


1.38 


10.0 


1.0 


1.23 


20.0 


2.0 


1.16 


30.0 


3.0 


1.12 


40.0 


4.0 


1.08 


50.0 


5.0 


1.10 


60.0 


6.0 


1.08 


70.0 


7.0 


1.07 


80.0 


8.0 


1.05 


90.0 


9.0 



RELATIVE SPEED 


ALFA 


UPLEG 


DOWNLEG 




8.5 


18.1 


0.6981 


15.6 


27.5 


0.6981 


24.8 


37.2 


0.6981 


34.4 


47.1 


0.6981 


44.2 


56.9 


0.6981 


54.1 


66.9 


0.6981 


64.0 


76.8 


0.6981 


74.0 


86.8 


0.6981 


83.9 


96.7 


0.6981 


7.7 


18.5 


0.7854 


14.7 


28.0 


0.7854 


24.0 


37.7 


0.7854 


33.7 


47.6 


0.7854 


43.5 


57.5 


0.7854 


53.4 


67.4 


0.7854 


63.3 


77.4 


0.7854 


73.3 


87.4 


0.7854 


83.2 


97.3 


0.7854 


6.8 


18.8 


0.8727 


13.9 


28.4 


0.8727 


23.2 


38.2 


0.8727 


33.0 


48.1 


0.8727 


42.8 


58.0 


0.8727 


52.7 


68.0 


0.8727 


62.7 


77.9 


0.8727 


72.6 


87.9 


0.8727 


82.6 


97.9 


0.8727 


6.0 


19.1 


0.9599 


13.1 


28.8 


0.9599 


22.6 


38.6 


0.9599 


32.3 


48.5 


0.9599 


42.2 


58.5 


0.9599 


52.1 


68.4 


0.9599 


62.1 


78.4 


0.9599 


72.0 


88.4 


0.9599 


82.0 


98.4 


0.9599 


5.2 


19.3 


1.0472 


12.4 


29.1 


1.0472 


21.9 


39.0 


1.0472 


31.7 


48.9 


1.0472 


41.6 


58.9 


1.0472 


51.6 


68.8 


1.0472 


61.5 


78.7 


1.0472 


71.5 


88.8 


1.0472 


81.5 


98.8 


1.0472 
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PROB . 

DETECTION 

1.0000 

0.9350 

0.8720 

0.8090 

0.7590 

0.7560 

0.7080 

0.7080 

0.6880 

1.0000 

0.9410 

0.8740 

0.8550 

0.8310 

0.7850 

0.8060 

0.7680 

0.7560 

1.0000 

0.9570 

0.9260 

0.9130 

0.9140 

0.8940 

0.8830 

0.8650 

0.8410 

1.0000 

1.0000 

0.9740 

0.9620 

0.9520 

0.9470 

0.9420 

0.9430 

0.9410 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 



r G. NO. 


SHIP 


SPEED 


RELATIVE SPEED 


ALFA 


PROB. 


LEGS 


SPEED 


RATIO 


UPLEG 


DOWNLEG 




DETECTI 


1.44 


10.0 


1.0 


4.3 


19.5 


1.1345 


1.0000 


1.24 


20.0 


2.0 


11.7 


29.4 


1.1345 


1.0000 


1.15 


30.0 


3.0 


21.4 


39.3 


1.1345 


1.0000 


1.12 


40.0 


4.0 


31.2 


49.2 


1.1345 


1.0000 


1.09 


50.0 


5.0 


41.2 


59.2 


1.1345 


1.0000 


1.08 


60.0 


6.0 


51.1 


69.2 


1.1345 


1.0000 


1.07 


70.0 


7.0 


61.1 


79.2 


1.1345 


1.0000 


1.05 


80.0 


8.0 


71.1 


89.2 


1.1345 


1.0000 


1.04 


90.0 


9.0 


81.0 


99.2 


1.1345 


1.0000 


1.55 


10.0 


1.0 


3.5 


19.7 


1.2217 


1.0000 


1.27 


20.0 


2.0 


11.1 


29.6 


1.2217 


1.0000 


1.18 


30.0 


3.0 


20.9 


39.5 


1.2217 


1.0000 


1.15 


40.0 


4.0 


30.8 


49.5 


1.2217 


1.0000 


1.11 


50.0 


5.0 


40.7 


59.5 


1.2217 


1.0000 


1.11 


60.0 


6.0 


50.7 


69.5 


1.2217 


1.0000 


1.08 


70.0 


7.0 


60.7 


79.5 


1.2217 


1.0000 


1.07 


80.0 


8.0 


70.7 


89.5 


1.2217 


1.0000 


1.07 


90.0 


9.0 


80.7 


99.5 


1.2217 


1.0000 


1.76 


10.0 


1.0 


2.6 


19.8 


1.3090 


1.0000 


1.36 


20.0 


2.0 


10.7 


29.8 


1.3090 


1.0000 


1.25 


30.0 


3.0 


20.5 


39.7 


1.3090 


1.0000 


1.18 


40.0 


4.0 


30.5 


49.7 


1.3090 


1.0000 


1.14 


50.0 


5.0 


40.4 


59.7 


1.3090 


1.0000 


1.12 


60.0 


6.0 


50.4 


69.7 


1.3090 


1.0000 


1.12 


70.0 


7.0 


60.4 


79.7 


1.3090 


1.0000 


1.08 


80.0 


8.0 


70.4 


89.7 


1.3090 


1.0000 


1.08 


90.0 


9.0 


80.4 


99.7 


1.3090 


1.0000 


2.19 


10.0 


1.0 


1.7 


19.9 


1.3963 


1.0000 


1.63 


20.0 


2.0 


10.3 


29.9 


1.3963 


1.0000 


1.44 


30.0 


3.0 


20.2 


39.9 


1.3963 


1.0000 


1.28 


40.0 


4.0 


30.2 


49.9 


1.3963 


1.0000 


1.27 


50.0 


5.0 


40.2 


59.9 


1.3963 


1.0000 


1.18 


60.0 


6.0 


50.2 


69.9 


1.3963 


1.0000 


1.15 


70.0 


7.0 


60.2 


79.9 


1.3963 


1.0000 


1.15 


80.0 


8.0 


70.2 


89.9 


1.3693 


1.0000 


1.15 


90.0 


9.0 


80.2 


99.9 


1.3963 


1.0000 
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MONTE CARLO TECHNIQUE 
PROBABILITY OF RANDOM ENCOUNTER 
NORMALIZED CONVOY ROUTE 



DEFINITE RANGE LAW 

DETECTION RANGE: 30.0 MILES 

LENGTH OF SUB PATROL LINE: 120.0 MILES 

X-DISTANCE TO PATROL LINE: 30.0 MILES 

SUBMARINE SPEED: 10.0 KNOTS 



NO. OF 


SAMPLES 


TAKEN 


FOR EACH 


SCENARIO 


>: 1000 




AVG. NO. 


SHIP 


SPEED 


RELATIVE SPEED 


ALFA 


PROB. 


LEGS 


SPEED 


RATIO 


UPLEG 


DOWNLEG 




DETECTION 


1.39 


10.0 


1.0 


14.1 


14.1 


0.0 


0.9212 


1.27 


20.0 


2.0 


22.4 


22.4 


0.0 


0.6714 


1.20 


30.0 


3.0 


31.6 


31.6 


0.0 


0.5720 


1.19 


40.0 


4.0 


41.2 


41.2 


0.0 


0.5512 


1.16 


50.0 


5.0 


51.0 


51.0 


0.0 


0.5324 


1.16 


60.0 


6.0 


60.8 


60.8 


0.0 


0.5221 


1.12 


70.0 


7.0 


70.7 


70.7 


0.0 


0.5000 


1.13 


80.0 


8.0 


80.6 


80.6 


0.0 


0.5181 


1.12 


90.0 


9.0 


90.6 


90.6 


0.0 


0.5133 
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PROBABILITY OF RANDOM ENCOUNTER 
EXTENDED MODEL 



DEFINITE RANGE LAW 

DETECTION RANGE: 30.0 MILES 

LENGTH OF SUB PATROL LINE: 90.0 MILES 

X-DISTANCE TO PATROL LINE: 30.0 MILES 

CROSSING ANGLE: 0.0 DEGREES (PERPENDICULAR) 

SUBMARINE SPEED: 10.0 KNOTS 



NO. OF 


SAMPLES 


TAKEN 


FOR EACH 


SCENARIO: 720 




AVG. NO. 


SHIP 


SPEED 


RELATIVE SPEED 


CROSS 


PROBABILITY 


LEGS 


SPEED 


RATIO 


UP LEG 


DOWNLEG 


POSIT 


OF DETECTION 


1.36 


10.0 


1.0 


14.1 


14.1 


45.0 


1.0000 


1.32 


15.0 


1.5 


18.0 


18.0 


45.0 


1.0000 


1.29 


20.0 


2.0 


22.4 


22.4 


45.0 


1.0000 


1.27 


25.0 


2.5 


26.9 


26.9 


45.0 


0.9861 


1.26 


30.0 


3.0 


31.6 


31.6 


45.0 


0.9250 


1.24 


40.0 


4.0 


41.2 


41.2 


45.0 


0.8528 


1.22 


50.0 


5.0 


51.0 


51.0 


45.0 


0.8139 


2.41 


90.0 


9.0 


90.6 


90.6 


45.0 


0.7250 


1.36 


10.0 


1.0 


14.1 


14.1 


50.0 


1.0000 


1.32 


15.0 


1.5 


18.0 


18.0 


50.0 


1.0000 


1.29 


20.0 


2.0 


22.4 


22.4 


50.0 


0.9847 


1.27 


25.0 


2.5 


26.9 


26.9 


50.0 


0.9375 


1.26 


30.0 


3.0 


31.6 


31.6 


50.0 


0.9069 


1.23 


40.0 


4.0 


41.2 


41.2 


50.0 


0.8444 


1.21 


50.0 


5.0 


51.0 


51.0 


50.0 


0.7958 


2.37 


90.0 


9.0 


90.6 


90.6 


50.0 


0.7056 


1.36 


10.0 


1.0 


14.1 


14.1 


55.0 


1.0000 


1.32 


15.0 


1.5 


18.0 


18.0 


55.0 


1.0000 


1.29 


20.0 


2.0 


22.4 


22.4 


55.0 


0.9292 


1.27 


25.0 


2.5 


26.9 


26.9 


55.0 


0.8792 


1.24 


30.0 


3.0 


31.6 


31.6 


55.0 


0.8333 


1.20 


40.0 


4.0 


41.2 


41.2 


55.0 


0.7792 


1.18 


50.0 


5.0 


51.0 


51.0 


55.0 


0.7500 


2.32 


90.0 


9.0 


90.6 


90.6 


55.0 


0.7056 
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AVG . NO . 


SHIP 


SPEED 


RELATIVE SPEED 


CROSS 


PROBABILITY 


LEGS 


SPEED 


RATIO 


UP LEG 


DOWNLEG 


POSIT 


OF DETECTION 


1.36 


10.0 


1.0 


14.1 


14.1 


60.0 


1.0000 


1. 32 


15.0 


1.5 


18.0 


18.0 


60.0 


0.9569 


1.28 


20.0 


2.0 


22.4 


22.4 


60.0 


0.8667 


1.24 


25.0 


2.5 


26.9 


26.9 


60.0 


0.7958 


1.21 


30.0 


3.0 


31.6 


31.6 


60.0 


0.7500 


1.17 


40.0 


4.0 


41.2 


41.2 


60.0 


0.6958 


1.15 


50.0 


5.0 


51.0 


51.0 


60.0 


0.6736 


2.26 


90.0 


9.0 


90.6 


90.6 


60.0 


0.6681 


1.36 


10.0 


1.0 


14.1 


14.1 


65.0 


1.0000 


1.32 


15.0 


1.5 


18.0 


18.0 


65.0 


0.9014 


1.26 


20.0 


2.0 


22 . 4 


22.4 


65.0 


0.7833 


1.22 


25.0 


2.5 


26.9 


26.9 


65.0 


0.7125 


1.19 


30.0 


3.0 


31.6 


31.6 


65.0 


0.6667 


1.16 


40.0 


4.0 


41.2 


41.2 


65.0 


0.6208 


1.14 


50.0 


5.0 


51.0 


51.0 


65.0 


0.6181 


2.23 


90.0 


9.0 


90.6 


90.6 


65.0 


0.6125 


1.36 


10.0 


1.0 


14.1 


14.1 


70.0 


1.0000 


1.31 


15.0 


1.5 


18.0 


18.0 


70.0 


0.8250 


1.24 


20.0 


2.0 


22.4 


22.4 


70.0 


0.7000 


1.21 


25.0 


2.5 


26.9 


26.9 


70.0 


0.6292 


1.18 


30.0 


3.0 


31.6 


31.6 


70.0 


0.5833 


1.15 


40.0 


4.0 


41.2 


41.2 


70.0 


0.5653 


1.13 


50.0 


5.0 


51.0 


51.0 


70.0 


0.5625 


2.23 


90.0 


9.0 


90.6 


90.6 


70.0 


0.5569 
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APPENDIX B 



MAIN COMPUTER PROGRAM FLOW DIAGRAM 
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yes (detection occurs) 




no 




yes 
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BASIC MODEL 




81 



EXTENDED MODEL 




no 



€ 



i 



SHIP SPEED <.25 KNO 1 



yes 



V 

INCREMENT 
SHIP SPEED 
BY 5 KNOTS 
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EXTENDED MODEL 




83 



SUBROUTINE WYEB 




84 



SUBROUTINE CPACHK 
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APPENDIX C 



************************* 

* NUMERICAL INTEGRATION * 

* BASIC MODEL * 

************************* 



CC^MCN’/ONE/SGN , YB , ZN 
INTEGER DALF 

DIMENSION UV(2)»W1V(2)»W2V(2) , V I V ( 2 I , V2V ( 2 ) , R IV ( 2 ) , 
1RNV(2),UW1V(2) ,UW2V(2) , TEMV ( 2 ) , WV ( 2 ) , VV ( 2 ) ,UWV ( 2 ) 

PI=3.1 4 15S26536 

FORMAT CONTROL COUNTER (INITIALIZE) 

I Z = 2 

SUB SPEED 

V = 10. 

SAMPLE SIZE 
JKL=960 

NP. CF ANGLE (ALFA) INCREMENTS 

I A 1=9 
NA=5 

NO. CF RHC INCREMENTS 
I BI=12 

DEFINE BOTH SUB VECTORS 
CALL VINS( V1V,0. ,V) 

CALL VINS( V2V,0. ,-V) 

SIZE OF DETECTION RADIUS (CEFINITE RANGE LAW) 

R=3C. 

LENGTH OF SUB PATRCL 
PI =1 20. 

INITIAL X-POSI TICN OF PATRCL LINE 

XP = R 

PRINT FIRST PAGE HEADING 
WPITE(6,60C)R,PL,XP,V, JKL 

INITIAL ANGLE 
ALF=0. 

LOOP TC INCREMENT ALFA 
DC 8C0 I A=1 , 1 A I 

CONVERSION OF RACIANS TO DEGREES 
I F ( I A — 1 ) 105,23,21 
23 DALF=ALF 
GC TC 22 

21 DALF= ( I A+3 ) *5 
INITIAL SHIP SPEED 

22 U=10. 

CALCULATE POSITION OF BCTTCM CF PATRCL LINE TO 
INSURE CROSSING AT CENTER OF PATROL 
YP=XP*TAN( ALF )-(PL/2. > 

C LOOP TO INCREMENT RhO 

DO 801 I B= 1 , I B I 
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C INITIALIZE LEG CCUNT 

SLEG=C. 

CCSA=CCS ( ALF ) 

SINA=SIN( ALF ) 

C DEFINE VECTOR OF SHIP VELOCITY 

CALL VINS(UVtU*COSAfU*SINA) 

C COMPLTE 1/V 

VI=1./V 

C COMPUTE SHIP/SUB SPEED RATIO 

RHC=U*VI 

C INITIALIZE SUB DIRECTION FOR DETERMINISTIC 

C CALCULATION 

Z N= 1 • 

INITIALIZE SUB POSITION FOR DETERMINISTIC 
CALCULATION 

YB = YP 

INITIALIZE DETECTION CCUNT 
D E T = C • 

COMPUTE MAXIMUM TIME FOR CALCULATIONS 
TMAX=( XP+P ) /(U*COSA J+l. 

DETERMINE RELATIVE SPEED VECTORS AND VALUES 
CF MAGNITUDES ANC ANGLES 
CALL VADD(V1V,UV»W1V,-1) 

CALL VADD(V2VtUV»W2Vt“l) 

CALL VVALIWIVtM *UW1V*BWI»111I 
CALL VVAL(W2V,W2,UW2V,BW2,111) 

DETERMINE IF DETECTION OCCURS 

EACH CCUNT GF THE LOOP INDEX REPRESENTS 
ONE SAMPLE 

DC 8 J=1 » JKL 

DETERMINE SUB POSIT 
INPUTS 

XP: X-D I STANCE TO PATROL LINE 

YP: Y-CCCRDINATE OF LOWER LIMIT CF PATROL 

LINE 

PL: LENGTH OF PATROL LINE 

OUTPUTS 

PC: RANGE TO SUB'S INITIAL POSIT 

BO: BEARING OF SUB'S INITIAL POSIT 

YP (VIA CCMMCN/ONE/I : INITIAL Y-CCORCINATE 

OF THE SUBMARINE 

R1V: INITIAL SUB POSITION VECTOR 
CALL WYEB< XP,YP,PL,RC,BC,R1V) 

C TIME TO ENC CF FIRST ( RANDCM ) LEG 







IF< SON) 700,105, 


702 














c 




TIME TO GC 


FOR 


SUB 


TO 


REACH 


LIMIT 


OF PATROL 


LINE 


c 




IF INITIAL 


SUB 


VELOCITY IS 


"DOWN" 








7 CC 


T 1= ( YB- YP ) *V I 

GO TO 703 
















c 




TIME TO GC 


FOP 


SUB 


TC 


REACH 


UPPER 


LIMIT CF 


PATROL 


c 




LINE 


















7 C 2 


T1=(YP+PL-YB)*V 


I 















7C3 T=T I 
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c 



INITIALIZE TOTAL TIME 



TN=0. 

C INDICATE THAT THE FIRST (RANDOM) LEG IS TO BE EX- 

C AMINE D 

MM=I 

C COMPUTE LENGTH OF RANDOM LEG 

GO TO 1 04 

C LOOP TO INCREMENT LEG CCUNT 

1 C 3 DC 1C IN=1,MM 

C CHECK TO DETERMINE INITIAL DIRECTION CF SUB 

IF(SGN)1C1 ,105,102 

C INITIAL DIRECTION IS UP; USE PARAMETERS FOR THE 

C UP LEG 

1C1 TW=TK2 

CALL COPI (LWV,UW2V) 

W=W2 

CALL COPI (WV,W2V) 

CALL COP I ( VV » V2 V ) 

GC TC 106 

INITIAL DIRECTION IS DOWN; USE PARAMETERS FOP THE 
DOWN LEG 

102 CALL C0PI(UWV,UW1V) 

TK=TW1 
W=W1 

CALL C0PI(WV,W1V) 

CALL C0PI(VV,V1V) 

REVERSE DIRECTION CF NEXT LEG 
106 SGN=-SGN 

COMPUTE RELATIVE LENGTH OF THIS LEG 
CALL VMULT (UWV»TEMV»TW) 

DETERMINE VECTOR FOR FINAL POSITION ON THIS LEG 
CALL VADD(R1V. TEMV,RNV,l) 

CALL CPACHK(R1V,RNV,UWV,RMIN, IN, J) 

INCREMENT LEG CCLNT 
SLEG=SLEG+1. 

DETERMINE IF DETECTION GCCURRED ON THIS LEG 
I F ( RM I N-P ) 3 , 3 , 1 1 

NO DETECTION THIS LEG; DETERMINE IF SHIP IS OUT 
OF RANGE BEYCND PATROL LINE 
11 TN=TN+T 
A IF(TN-TMAX)6,8 ,8 

IF THIS IS NOT THE FIRST (RANDOM) LEG, CONTINUE 
6 IF(MM-1)105,9,10 
10 CONTINUE 
GO TO 8 

LENGTH CF TIME REQUIRED FOR SUB TO TRAVEL ENTIRE 
SPAN OF PATROL LINE 
9 T=PL*VI 

C PROVIDE FOR POSSIBILITY OF 10000 LEGS MAXIMUM 

MM=1 0000 

C COMPUTE LENGTHS CF SECOND AND SUCCEEDING LEGS 

104 TW1=T*W1 
TW2=T*W2 
GC TC 103 
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c INCREMENT CETECTICN COUNTER 

3 CET=DET+1 • 

8 CONTINUE 
XJKL=JKL 

C COMPUTE PROBABILITY OF DETECTION 

P=DET /X JKL 

C COMPUTE AVERAGE NO. OF LEGS FOR THIS SCENARIO 

XLEG=SLFG/XJKL 

C PRINT RESULTS FOP EACH SPEED RATIO 

WRITE !6 ,6C 1 ) XLEG,U,RHC, WI ,W2,DALF,P 

C INCREMENT SHIP SPEEC BY 5 KNOTS IF RHC<7, 

C THEN INCREMENT BY 10 KNOTS 

I F ( 16-7)16,17,17 

16 U=U+5. 

GC TC 801 

17 U=U+ 10* 

601 CONTINUE 

C FORMAT OUTPUT TO INCLUDE RESULTS FOR TWO 

C VALUES OF ALFA ON FIRST PAGE, AND THREE 

C VALUES OF ALFA CN SUBSEQUENT PAGES 

IF! IA-IZ)12, 13,105 

C ONE LINE FEEC BETWEEN VALUES OF ALFA 

12 WR I TE ( 6 ,602 ) 

GC TO 20 

C BEGIN NEW PAGE; INSERT COLUMN HEADINGS 

13 I Z=I Z + 3 

I F ( IA-9)19,800,105 

19 WP I TF ( 6 , 603 ) 

C INCREMENT ALFA BY 25 DEGREES CN THE FIRST 

C ITERATION, THEN BY 5 DEGREES ON SUBSEQUENT 

C iterations 

20 ALF=ALF+NA*(PI/36.) 

800 NA=1 



6CC FCRMAT11H1.15!/) ,T16, 'PROBABILITY OF RANDOM ENCOUNTER' 
1/T26, 'BASIC MODEL '///T6, 'DEFINITE RANGE L AW • /T8 , • DETE • 

2, 'CTI0N RANGE: ',F5.1,' MILES '//T6, • LENGTH OF SUB PAT* 

3, 'ROL LINE: *,F5.1,' MILES' //T6 ,* X-DISTANCE TO PATROL* 
5,' LINE: ' , F5. 1 , ' MILES' //T6, 'SUBMARINE SPEEC :', F5. 1, 
5' KNOT S'//T6,'NC. OF SAMPLES TAKEN FOR EACH SCENARIO:' 
fc,' • ,I4///1X, • AVG.NO. SHIP SPEED RELATIVE SPEEDS • 

7, 'ROUTE PROBABILITY 0F'/2X,'LEGS SPEED RATIO UPLE' 

8, 'G OOWNLEG ANGLE DETECTION'//) 

601 FORMAT! 1X,F5.2»1X,2!2X,F4.1 ) , 2 < 3X , F5. 1 ) , 6X , I 2 , 8X, F6.4 ) 

602 FORMAT!/) 

603 FCRMAT(1H1,15(/),1X,'AVG.N0. SHIP SPEEC RELATIVE SP' 

1, 'EEDS ROUTE PROBABILITY OF • /2X, 'LEGS SPEED RATI' 

2, 'C UPLEG DCWNLEG ANGLE DETECTION'//) 

105 STOP 
END 
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************************* 

* NUMERICAL INTEGRATION * 

* EXTENDED MODEL * 

************************* 



C 

C 

c 

c 



COMMON /ONE /SGN , Y8 » ZN 

DIMENSION UV(2),W1V(2),W2V(2)*V1V(2),V2V(2),R1V(2), 
1RNVI2) ,UW1V<2) ,UW2V(2) ,TEMV ( 2 ) , WV < 2 ) , VV < 2 ) , UWV< 2 1 

PI=3.1415926536 

FORMAT CONTROL COUNTER UNITIALIZE) 

I Z=2 



V=10. 



SUB SPEED 



SAMPLE SIZE 
JKL=960 

NO. CF TRACK CROSSING INCREMENTS 
I AI=10 

NO. OF RHO INCREMENTS 
I B I = 1 1 

DEFINE BOTH SUB VECTORS 
CALL VINSt V1V,0. ,V) 

CALL VINSI V2V,0. .-V ) 

SIZE OF DETECTION RADIUS (CEFINITE RANGE LAW) 

R=30. 



LENGTH OF SUB PATROL 
PL=120. 

INITIAL X-POSITICN OF PATRCL LINE 

XP=R 



PRINT FIRST PAGE HEADING 
WRITE(6,600)R,PL,XP,V,JKL 

INITIAL ANGLE 
ALF=0. 

YP=-PL/2. 

LOOP TO INCREMENT CROSSING POSITION 
DO 800 I A=1 « I AI 



U=10. 



INITIAL SHIP SPEED 



CALCULATE POSITION OF BOTTCM OF PATRCL LINE TO 
INSURE CROSSING AT CENTER OF PATROL 

LOOP TO INCREMENT RHO 
DO 801 I B = 1 y I B I 

DELETE DATA FOR 5<RH0<9 
IF(U-50. 120,20,21 

GENERATE DATA FOR RH0*9 
21 IFIU-90. >17,22,105 

INITIALIZE LEG COUNT 
20 SLEG*0. 
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22 CCSA=CCS( ALF ) 

S INA=SI N( ALF ) 

C CEFINE VECTOR OF SHIP VELOCITY 

CALL VINS(UV,U*COSA ,U*SINA) 

C COMPUTE 1/V 

VI=i./V 

C COMPUTE SHIP/SUB SPEEC RATIO 

RHO=U*V I 

C INITIALIZE SLB DIRECTION FOR DETERMINISTIC 

C CALCULATION 

ZN = 1. 

C INITIALIZE SUB POSITION FOR DETERMINISTIC 

C CALCULATION 

YP=YP 



C 

C 



c 

c 



c 

c 



c 

c 



c 

c 



INITIALIZE DETECTION COUNT 

C E T = 0 • 

COMPUTE MAXIMUM TIME FOR CALCULATIONS 
TMAX=(XP+R)/(U*COSA)+l. 



AND VALUES 



DETERMINE RELATIVE SPEED VECTORS 

OF MAGNITUDES ANC ANGLES 
CALL VADD(VlV»UV»WIVt-l) 

CALL VADD(V2V»UV*V*2V»-1 ) 

CALL VVAL(V.IV,W1,UW1V,BW1,111) 

CALL VVAL(W2V,W2,UW2V,BW2,111) 

DETERMINE IF DETECTION OCCURS 

EACH COUNT OF INCEX J REPRESENTS ONE SAMPLE 
DO 8 J=I » JKL 

DETERMINE SUB POSITION 
INPUTS 

XP: X-D I STANCE TO PATROL LINE 

YP: Y-CCGRDINATE OF LOWER LIMIT OF PATROL 

LINE 

PL: LENGTH OF PATROL LINE 

OUTPUTS 

RC: RANGE TO SUB'S INITIAL POSITION 

BO: BEARING OF SUB'S INITIAL POSITION 

YB (VIA CCMMCN/ONE/) : INITIAL Y-COCRDINATE 

OF SUBMARINE 

PIV: INITIAL SUB POSITION VECTOR 
CALL WYEB(XP»YP»PLfRCfBC»RlV) 

TIME TO END OF FIRST (RANDOM) LEG 
IF(SGN)70C»105»702 

TIME TO GO FOR SUB TO REACH LIMIT OF PATROL LINE 
IF INITIAL SUB VELOCITY IS "DOWN" 



7CC Tl=< YB-YP)*VI 
GO TO 703 



TIME TO GO FOR SUB TO REACH UPPER LIMIT OF 
LINE 



PATROL 



7C2 Tl=( YP+PL-YB)*VI 
703 T=T I 



TMO. 



INITIALIZE TOTAL TIME 



C 

C 



MM=1 



INDICATE THAT THE FIRST (RANDOM) LEG IS TO BE EX- 
AMINED 
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c 

c 



103 



101 



102 



106 



11 

4 



6 

10 



C 

C 



104 



COMPUTE LENGTH OF RANDOM LEG 
GO TO 104 

LOOP TO INCREMENT LEG COUNT 
DO 10 IN=1,MM 

CHECK TO DETERMINE INITIAL DIRECTION OF SUB 
IF(SGN)101, 105,102 

INITIAL DIRECTION IS UP; USE PARAMETERS FOR THE 
UP LEG 

TM=TW2 

CALL COPI(UWV,Uk2V) 

W S W2 

CALL COP I ( WV » W2V ) 

CALL COP I C VV* V2V ) 

GO TO 106 

INITIAL DIRECTION IS DOWN; USE PARAMETERS FOR THE 
DOWN LEG 

CALL C0PI(UWV,UW1V) 

TW=TW1 

W=W1 

CALL C0PI(WV,W1V» 

CALL C0PI(VV,V1V) 

REVERSE DIRECTION CF NEXT LEG 
SGN=-SG N 

COMPUTE RELATIVE LENGTH OF THIS LEG 
CALL VMULT(UWV,TEMV,TW! 

DETERMINE VECTOR FOR FINAL POSITION ON THIS LEG 
CALL VADD(R1VjTEMV»RNV» 1 ) 

CALL CPACHK(RlV*RNV»UWV»RMIN»INtJ) 

INCREMENT LEG CCUNT 
$LEG*SLEG+1* 

DETERMINE IF DETECTION OCCURRED ON THIS LEG 
IF(RMIN-RJ3,3,11 

NO DETECTION THIS LEG; DETERMINE IF SHIP IS OUT 
OF RANGE BEYOND PATROL LINE 
TN=TN+T 

IF(TN-TMAX)6»8*0 

IF THIS IS NOT THE FIRST (RANDOM! LEG* CONTINUE 
I F( MM— 1 >105*9*10 
CONTINUE 
GO TO 8 

LENGTH OF TIME REQUIREO FOR SUB TO TRAVEL ENTIRE 
SPAN OF PATROL LINE 
T»PL*VI 

PROVIDE FOR POSSIBILITY OF 10000 LEGS MAXIMUM 
MM=*10000 

COMPUTE LENGTHS OF SECOND AND SUCCEECING LEGS 
TW1=T*W1 
TW2*T*W2 
GO TO 103 

INCREMENT DETECTION COUNTER 
DET*DET+1* 

CONTINUE 

XJKL^JKL 
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c 

c 



c 

c 



c 

c 



COMPUTE PROBABILITY CF DETECTION 
P=DET/XJKl 

COMPUTE AVERAGE NO. CF LEGS FOR THIS SCENARIO 
XLEG=SLEG/XJKL 

GENERATE CROSSING COORDINATE FCR PRINTOUT 
X YP=— YP 

PRINT RESULTS FOR EACH SPEED RATIO 
WRITE(6»601) XLEG»U»RHC»W1 ,W2,XYP,P 

INCREMENT SHIP SPEEC BY 5 KNCTS IF RHCO, 

THEN BY 10 KNCTS 
IF< 16-5)16,17,17 

16 U=U+5. 

GC TC 801 

17 U=U+1C. 

8C1 CONTINUE 

FORMAT OUTPUT TO INCLUDE RESULTS FCR TWC 
VALUES OF ALFA ON EACH PAGE 
IFIIA-IZH2, 13,105 

ONE LINE F EEC BETWEEN VALUES OF ALFA 
12 WRITE (6 ,602) 

GO TC 18 



BEGIN NEW PAGE; 
13 I Z= I Z+2 

I F < I A— 10)19,18,105 
19 WPI TE ( 6 ,603 I 



INSERT COLUMN HEADINGS 



INCREMENT CROSSING POSITION BY 5 MILES 
UNTIL WITHIN 30 MILES OF PATROL TERMINUS, 

THEN INCREMENT BY 20 MILES 
19 I F ( I A— 7)14,15,15 

14 YP=YP-5. 

GO TO 800 

15 YP=YP-20. 

800 CONTINUE 

6CC FORMAT ( 1H1 , 8 ( / ) , T16 , 'PROBABILITY CF RANDOM ENCOUNTER* 
1/T2A, 'EXTENDED MODEL* 

1///T6, 'DEFINITE RANGE L AW • /TB ,* DETECT ION RANGE: *.F5.1 

2, ' MILES'//T6, 'LENGTH CF SUB PATROL LINE: ',F5.1,* MI* 

3, 'LES'//T6,'X— DISTANCE TO PATROL LINE: ',F5.1,' MILES' 
A//T6, 'CROSSING ANGLE: 0.0 DEGREES ( PERPENC ICULAR ) • 
5//T6, 'SUBMARINE SPEED :• »F5. 1, • KNCTS • //T6 ,' NO • OF SAM* 

6, 'PLES TAKEN FOR EACH SCENARIO: • , 14/ // 1 X , • AVG.NO. SH' 

7, 'IP SPEED RELATIVE SPEEDS CROSS PROBABILITY OF'/ 
82X, • LEGS SPEED RATIC UPLEG DOWNLEG POSIT DE' 
9 , • TECTI ON • // ) 

601 FCRMAT(1X,F5.2,1X,2(2X,FA.1 ),2(3X,F5.1),5X,F5.1,7X, 
1F6.4 ) 

602 FORMAT!/) 

503 F0RMAT(1H1,17( /) .lX.'AVG.NO. SHIP SPEED RELATIVE SP' 

1, 'EEDS CROSS PROBABILITY 0F'/2X,'LEGS SPEEO RATI* 

2, 'C UPLEG DOWNLEG POSIT DETECTION'//) 

105 STOP 
END 
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******************************** 

* NUMERICAL INTEGRATION MODELS * 

* SUBROUTINES * 

******************************** 



SUBROUTINE WYEB< XP,YP ,PL,RC, BC,A) 

THIS SUBROUTINE GENERATES SECUENTI AL Y- 
COORD I NATES ON A SPECIFIED LINE SEGMENT; 
WHEN AN ENC PCINT OF THE LINE IS GENERATED, 
THE DIRECTION GF THE SEQUENCE IS REVERSED 

INPUTS 

Y6=Y-P0S I T ION CF SUB ON PATROL LINE 
SGN= INITIAL DIRECTION OF SUB 



CCMMON/ONE/SGN,YB,ZN 
DIMENSION A ( 2 ) ,B(2) 

INCREMENT SU8 POSITION BY 500 YDS IF ZN=1; 
DECREMENT SUB POSIT BY 500 YDS IF ZN=-I 
YB=YB40.25*ZN 

REVERSE SUB COURSE WHEN IT REACHES ENC CF 
PATROL LINE 

IF(ZN*(2.*(YP-Ye > + PL) + PU9,9, 10 

POSITION SUB EXACTLY AT TERMINUS OF PATROL 
LINE PRIOR TC TURN 
9 YB=YP+0*5*PL*( 1 • + ZN ) 

C REVERSE COURSE CF SUB 

ZN=-ZN 
10 SGN=ZN 

C CETERMINE INITIAL RANGE AND BEARING TC SUB 

A ( 1 ) =XP 
A ( 2 ) =Y5 

CALL VVALI A, R0,B, 80,101 ) 

RETURN 

END 
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Sl'BRCUT l NE VINS(A,X,Y) 



C 

C 



THIS SUBROUTINE INSERTS 
COMPONENTS INTC ARRAY A 



THE X ANC Y VECTOR 



DIMENSION 
A ( 1 ) = X 
A C 2 1 =Y 
RETURN 
END 



aid 



SUBROUTINE COPI(A,B) 



C 

c 



THIS SUBROUTINE GENERATES THE VECTOR A SUCH THAT 
FOR A GIVEN VECTOR B, A = B 



DIMENSION A ( 1 1 , B ( 1 ) 
DC 1 1*1.2 
A ( I I =B( I ) 

RF T UPN 

END 



SUBROUTINE VADD ( A , B , R , N ) 



THIS SUBROUTINE PERFORMS THE FOLLOWING VECTOR 
OPERATIONS: 

IF N>C : R = A + B (VECTOR ADCITICNI 
IF N<0: R = A-B (VECTOR SUBTRACT I ON ) 

IF N=0 : THE DOT PRODUCT OF A AND B IS 
RESULT IS PLACED IN THE FIRST MEMBER 



FORMED AND 
CF ARRAY R 



DIMENSION A( 1 1 ,6(1) ,R(1) 

IF(N)l,2,3 

ADD VECTORS 
3 DO 4 1*1,2 
A R(I)=A(n+B(I) 

RETURN 

SUBTRACT VECTCRS 

1 DC 5 1*1.2 

5 R(I)=A(n-e(I) 

RF T URN 

COMPUTE DCT PRODUCT CF VECTORS 

2 DO 6 1=1,2 

6 R(I)*A( I )*B< I > 

R ( 1 ) =P < 1 ) +P ( 2) 

RFTUPN 

END 
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SUBROUTINE VMULT(A,B,X) 



THIS SUBROUTINE MULTIPLIES THE VECTOR A BY THE 
SCALAR X TC FORM THE VECTOR B 

DIMENSION A (1 ) * B ( 1 ) 



SUBROUTINE VVAL< A , VL ,UV , ANGLE ,M ) 



THIS SUBROUTINE DETERMINES THE FOLLOWING 
VALUES OF At 

A=INPUT VECTOR (SUPPLIED BY THE CALLING PROGRAM ) 

VL*MAGNITUOE OF A 

UV=UNI T VECTOR PARALLEL TO A 

ANGLE =ANGLE OF A MEASURED CCW FROM X-AXIS 



p»n -l t — 1 o 




END 





RETURN 



C 




C 



COMPUTE BOTH LENGTH AND ANGLE OF VECTOR 

5 VL=SQRT(A(1)**24A<2»**2> 



GO TO 2 



C 



GENERATE UNIT VECTOR ONLY 
6 VL=$QRT ( A ( 1 1 **2+A ( 2 1**2 ) 



nn i t-i o 




c 



COMPUTE ANGLE AND GENERATE UNIT VECTOR 



7 ANGLE=ATAN2( A ( 2 ) , A(1J J 
GO TC 6 
END 
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SUBROUTINE C P ACHK I R 1 , RN , U WV , RCP A , I N , J ) 



THIS SUBROUTINE DETERMINES THE CLOSEST PCINT 
TO THE ORIGIN OF ANY LINE SEGMENT. THE ENO 
POINTS OF ViH I CH ARE DEFINED BY VECTORS 

R1=RFFEPENCE VECTOR 

RN=TERMINAL VECTOR 

UW V=UN I T VECTOR PARALLEL TO CRM 

RCPA=DI STANCE AT CLOSEST PCINT OF APPROACH (CPA) 
I N=LEG NO. BEING EXAMINED AT THIS EXECUTION 
J = S A MPLE NO. BEING EXAMINED AT THIS EXECUTION 



DIMENSION PHI) ,RN( 1 ) ,UWV(1 > ,X( 2),Y(2),CPA(2),Z(2), 

IP T ( 2 ) 

CCMMCN/TWC/Y.N 

C COMPUTE DCT PRCCLCT: R1*UWV=X 

CALL VADD(PI»UWV»X,0) 

IF(X( 1) ) 1,2,2 

C CPA WILL OCCUR ON THIS LEG; COMPUTE THE COT FROD- 

C LCT : RN*UWV=X 

1 CALL VADD(RN»UWV»X,0) 

I F ( X ( 1 ) )3,3,4 

C CPA WILL NOT CCCUR ON THIS LEG; CPA=R1 

2 CALL VVAL(P1,RCPA,X,D,100) 

RETURN 

C CPA WAS NOT REACHEO ON THIS LEG; CPA=RN 

3 CALL VVAL(PN,RCPA,X, 0,100) 

RETURN 

C DETERMINE CPA; ASSUME LEG HAS LENGTH < 250 MILES 

4 DO 5 1=1, 5C0 

C USE 2-MILE INCREMENTS 



P=2. * I 



C 



INCREMENT CRM BY 2 MILES 
CALL VMULT ( UWV » Y , P ) 



C 



ADD TC REFERENCE VECTOR 
CALL VADD (P1»Y»X,1) 



C 



COMPUTE DCT PRODUCT X*UWV=Z 
CALL VADD( X »UWV» Z ,0) 



C 

c 



I F ( Z ( 1) ) 5, 6 , 7 



IF Z BECOMES POSITIVE, MAKE FINE INCREMENT CHECK 
BELOW; IF Z(1)=0, X=CPA; OTHERWISE TRY AGAIN 



C 



INCREMENT CRM BY 0.1 MILE 



C 



PRESERVE PRECECING VECTOR 
7 CALL VMULT(UWV,Y,P-2. ) 

CALL VADC(P1 ,Y,X, 1) 



C 



CHECK LAST 2 MILES OF DRM IN 0.1-MILE INCREMENTS 
DC 8 1=1,21 
P=0. 1*1 



C 



INCREMENT CRM BY 0.1 MILE AND FIND CPA 
CALL V M ULT ( UWV , Y , P ) 

CALL VADD(X,Y,RT,1) 



IF LL VAUUl A § T I fll 

CALL V ADD ( RT ,UWV , Z , 0 ) 
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ERROR EXISTS IN PROGRAM IF NEXT INSTRUCTION IS 
EXECUTEO 

WRITE (6 ,61 ) I N, J 

61 F0RMATC//10X, 'LEG',1X,I4,1X,'CF SAMPLE' , 15, IX, • GENER' , 
l'ATED ERROR IN CPACHK ROUTINE'//) 

IFIZC 1) ) 5,6,7 
5 CONTINUE 

DRM LEG>250 MILES IF NEXT INSTRUCTION IS 



60 X,'CF SAMPLE', 15, IX, 'IS 




RETURN 



C 



COMPUTE CPA DISTANCE 



6 CALL VVAL(X,RCPA,Y,C, 100) 
RETURN 
END 
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APPENDIX D 



r **♦*♦*<[**** i* <;<:******* ** 

r * c A N P P M Sampling * 

* MIDPOINT CRHSSING * 

f *★ * * *** ****** * A ** **** ** 

cr mmon/ol r / s n n ♦ Y b 

DIMENSION I IV (2) »WlV(2)tW2V(2 I , V 1 V ( 2 ) , V2V ( 2 I ,R1V( 2 ) , 
]PNV( 2I.UW1 V(2I ,UW2V(2 ) ,TEMV( 2),WV<?) , VV( 2 I , UWV ( 2 ) 

C INITIH ! ?F URN GENFRATQP 

nr. SI T=HPN( C ) 

PT=3. 1 '-I c 9?^ t ?3^ 

r sup SPfFn 

V = I A . 

r S AMr If c I 7 F 

,)K( *ir v - 0 

C Mf. OF ANGLE (ALFA) INCREMENTS 

I 4 1 = 17 

0 NO. r>F F HO INC RFMFNTS 

T P I =0 

C DEFINE D OTH SUB V c CTCPS 

CALL VJMf ( VIV.O. »V) 

CALL VTNS( V2V,0. ,-V) 

r S I 7 r OF OETFCTION PAOILJS (DEFINITE RANGE LAW) 

P =3C. 

C LENGTH of sup PATPOL 

PL =1 20. 

C JMITIM X-PPSITION OF PAT^pl LTNE - 

xn=F 

WFITF(S,FCR)R,RL,XP,V t JKL 
I N I T I A L ANGl c 
A I F = C.O 

C l TP TO INCREMENT ALFA 

r r POO IA = 1 ,TAI 
g INITIAL SHIP SPEED 

• i = !0. 

CALCULATF POSITION OF POTTO* OF »ATR°L LINF TO 
INSURF CROSSING AT CENTER OF PATROL 

YF=¥F*TAN( ALF)-(PL/2. ) 
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C LOOP TO INCREMENT RHO 

DO 801 I B=1 t I B I 

C INITIALIZE LEG COUNT 

SLEG=0. 

C0SA=C0S( ALF) 

SI NA=SI N( ALF I 

C DEFINE VECTOR OF SHIP VELOCITY 

CALL VINS(UV f U*COSAtU*SINA) 

C COMPUTE SHIP/SUB SPEED RATIO 

RH0=U/V , 

INITIALIZE SUB DIRECTION FOR DETERMINISTIC 
CALCULATION : 

SGN=1. 

INITIALIZE DETECTION COUNT 



DET=0. 

COMPUTE MAXIMUM TIME FOR CALCULATIONS 
TMAX=(XP+R)/(U*C0SA)+1. 

DETERMINE RELATIVE SPEED VECTORS AND VALUES OF 
MAGNITUDES AND ANGLES 

CALL VADD(VlV f UV»WlVf-l) 

CALL VADD( V2V» UV» W2V »-l ) 

CALL VVAL(W1V,W1,UW1V,BW1,111) 

CALL VVAL ( W2Vf W2 » UW2V « BW2 till) 

DETERMINE IF DETECTION OCCURS 

EACH COUNT OF THE LOOP INDEX REPRESENTS ONE SAMPLE 
DO 8 J=1 * JKL 

DETERMINE SUB POSIT (UNIFORMLY DISTRIBUTED) 

INPUTS 

XPi X-D I STANCE TO PATROL LINE 
YP: Y-COORDINATE OF LOWER LIMIT OF PATROL 

LINE 

PL: LENGTH OF PATROL LINE 

OUTPUTS 

R0: RANGE TO SUB'S INITIAL POSIT 

BO: BEARING OF SUB'S INITIAL POSIT 

YB (VIA COMMON/ONE/): INITIAL Y-COORDINATE 

SUB 

RIV: INITIAL SUB POSITION VECTOR 
CALL WYEB(XP,YP,PL,R0,B0,R1V) 

TIME TO END OF FIRST (RANDOM) LEG 
IF(SGN)700,105,702 

TIME TO GO FOR SUB TO REACH LIMIT OF PATROL LINE 
IF INITIAL SUB VELOCITY IS "DOWN" 

700 Tl=( YB-YP ) /V 
GO TO 703 

TIME TO GO FOR SUB TO REACH UPPER LIMIT OF PATROL 



101 



r 



l INF 



7 r? t i = ( yp + fl -yr ) t v 

T rT' 

C INITIALIZE rnm TIMP 

Tr =c. 

r INDICATE that THE FIRST (RANDOM) LFG is TO of 

r ’■xahinet 

MM=1 

C COMPUTE length of random leg 

or to 104 

r LCCR T 0 INCREMENT LEG COUNT 

10’ or 1C IN=1,MM 

r CHECK TO DETERMINE initial direction of sur 

T F( SGNM ri ,105,102 

C INITIAL DIRECTION IS IJP; US C PARAMETERS C 0R THF 

C UP LFG 

I'M JW=TW? 

CALL COP I ( LIWV , LJW2 V) 

W= W2 

C/Ll COPI ( WV, W2V) 

CALL COP I ( VV,V2V) 

GC TP 106 

C INITIAL DIRECTION IS DOWN; USE PARAMETERS FOR THE 

C DOWM LFG 

102 CALL Cnpl I UWV,UW1 V) 

TH=TW1 
W = W 1 

CALL Cnpl(WV,WlV) 

CALL C0PI(VV,V1V) 

C PEVPFSE OIRFCTION OF NEXT L C G 

ICO SCN=-SGN 

C COMPUTF RELATIVE LENGTH O c THIS LEG 

CALL VMUL T ( IJVIV , TE M V , TW I 

C DETERMINE VECTOR FOR FINAL POSITION ON THIS LEG 

CALL VADD(P1V,TEMV,RNV,1> 

CALL CPACFK ( RlV,RNV,UWV,RMIN , IN , J ) 

TMCRFMFNT leg count 

SL EG = SLEG+ 1 . 

C DETPRHINF IF DETECTION OCCURRED ON THIS LEG 

IP (PM! N— P ) 3 , 3 , l t 
11 tlsTN+t 

•a !F( TN-TMAX 16.5 ,8 

6 IFLHM-i Jiff ,0,10 

1 1 CCNTINUF 
GC TO R 
i T = oi/V 
MM»1 COCO 
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104 TW1=T*W1 
TW2=T*W2 
GO TO 103 

3 DET=DET+1. 

8 CONTINUE 
X JKL= JKL 
P=DET/XJKL 
XLEG=SLEG/XJKL 

WRITE ( 6 , 601 ) XLEG,U,RHC,W1 ,W2,ALF,P 
801 U=U+10. 

WP ITE ( 6 , 602 ) 

602 FORMAT ( / ) 

800 ALF=ALF+PI/36. 

600 F0RMATUH1/////T16, 'PROBABILITY OF RANDOM ENCOUNTER'./ 
1//T6, 'DEFINITE RANGE LAW • /T8, 'DETECTI ON RANGE: '.F5.1. 
2' MILES'//T6, 'LENGTH OF SUB PATROL LINE: '.F5.1. 1 MIL 1 
3»'ES'//T6, 'X-DI STANCE TO PATROL LI NE: • , F5. 1, ' MILES' // 
4T6, 'SUBMARINE SPEED: *.F5. 1. • KNOTS • //T6, • NO. OF SAMPL » 

5, 'ES TAKEN FOR EACH SCENARIO: ',I4///T3,' AVG. NO. SH' 

6, 'IP SPEED ' 

7, 'RELATIVE SPEED ALFA PROB. DETECT ION • /T5 , 'LEGS ' 

8, ' SPEED RATIO UPLEG OOWNLEG'// ) 

601 FORMAT! T4, F5. 2 ,2X , 3 (1 X, F5.1 ) ,2X,F5. I , 3X, F6.4,7X, F6.4) 

105 STOP 
END 



C 

C 

c 

c 

c 



SUBROUTINE WYEBI XP , YP , PL ,R0, BO, A) 

INPUTS 

YB=Y-POSI TION OF SUB ON PATROL LINE 
YB IS UNIFORMLY DISTRIBUTED 

SGN=INITIAL DIRECTION I EQUALLY LIKELY) OF SUB 
VELOCITY 

COMMON/ONE/SGN, YB 
DIMENSION A ( 2 ) , B( 2 ) 

DETERMINE SUB POSIT RANDOMLY 

YB=YP+URN<1)*PL 

DETERMINE INITIAL DIRECTION 

A I F (URN( 1 )-0. 5) 5 ,6 ,7 
6 GO TO A 



5 SGN=1. 

GO TO 8 

DOWN 

7 SGN=-1 . 

DETERMINE INITIAL RANGE ANO BEARING TO SUB 

8 A ( 1 ) =XP 
A(2)*YB 

CALL VVAL( A, RO, B,80,101 ) 

RETURN 

END 
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St’RpruriNr CRACHK ( Rl, PN,tJWV,RCPA, IN, J) 

ri M rLSION nill),RW(ll ,UWV III f X(2),Yl2l,CPA(?l,ZI2l 

r r M v r n / t wo / v ,n 

r m=»f c frence vector 

r. p^‘=t c pminal vectgp 

C UK V = UN I T V E C T O R CF DRM 



r COMPIITI HOT PRODUCT: Rl*tJWV=X 

call vAnn(Ri,uwv,x,o) 
if(x( m i ,?, 2 

C CPA Wilt OCCUR ON THIS LEG: COMP'JTF the dot 

f PRODUCT: RN*UWV=X 

1 C/Ll VADD(PN,UWV,X,0) 

I p ( X ( 1 1 ) 3 , 3 ,4 

C CPA Wilt NOT OCCUR ON THIS LEG: CPA=P1 C™°n T F 

c CP A 

? call vva li pi ,rcpa ,x» o, loo ) 

Rf T UFN 

C CPA WAS NOT REACHED ON THIS LEG: CPA=RN: CPMPUTF 

C C r A 

3 CALt VVAl (FN,RCPA,X,D,100) 
return 

C DETERMINE CPA: ASSUME LEG HAS LFNGTHC 250 MILES 

i nr S 1=1, ECO 

r usf half-mile incrfmfnts 



r 

r 



p=l /2. 

INCPFMFMT DP M RY 1/2 mile 
C A t l VMUl T(UWV,Y,P) 

ADD TO PEFFRFNCE VECTOR 
CALL VADDIBl ,Y,X, 1) 

COMPIITF DOT PRODUCT X DOT IJWV=Z 
CALI VADrt X,UWV,Z ,0) 



IF 7 CHANGES SIGN, OR 
OTHERWISE TRY AGAIN 



R6C0MFS ZFRO VECTOR X=CPA: 



TF< Z(ll ) 
CONTINUE 



5,5, h 



PPM L EG>2 50 MILES 
WRITE ( ft , AC ) IN, J 

FOPMATt //1CX, 'LEG* , 1 X , 1 4 , l X , ‘OF SAMPL E' ,15, IX, 'IS 
1 GREATER THAN 250 MILES'//) 

RF TURN 

a CALL VVAl (X,RCPA,Y,D, 100) 

Pf TURN 
END 
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c 

c 



SUBROUTINE VINS(A,X,Y> 

THIS SUBROUTINE INSERTS THE X AND Y VECTOR 
COMPONENTS INTO ARRAY A 




END 



SUBROUTINE COPI(A,B) 

THIS SUBROUTINE GENERATES THE VECTOR A SUCH THAT 
FOR A GIVEN VECTOR B, A=B 



SUBROUTINE VADDI A , B, R , N) 

THIS SUBROUTINE DOES THE FOLLOWING VECTOR 
OPERATIONS: 

IF N>0 : R=A+B 
IF N<0: R=A- B 

IF N=0 : THE DOT PRODUCT OF A AND B IS FORMEO AND 
RESULT IS PLACED IN THE FIRST MEMBER OF ARRAY R 

DIMENSION Am f B<n,Rm 
IF(N)1»2»3 

C ADD VECTORS 

3 DO 4 1=1.2 

4 R(n=A(n+B(I) 



uu i 1=1. 
1 A( I )=B( I ) 
RETURN 
END 



DIMENSION All), Bill 
DO 1 1=1,2 



C 



1 UU D I-ili 

5 RII)=AU)-BII) 



RETURN 

SUBTRACT VECTORS 
1 DO 5 1=1,2 



RETURN 



C 




COMPUTE DOT PRODUCT OF VECTORS 
2 DO 6 1=1,2 
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SUPPriJTIMf VMULT(4,R,X) 



T H lc C(|ppdi)T I NE MULTIPLIES THE V EC T OP A BY the 
SCALAR X TO FORM THE VECTOR B 

niMF»'CirN A ( 1 ) » B ( 1 ) 

or 1 1=1.2 

P ( I ) = X*A ( I ) 

p r t upn 

r no 



c I T E CUT INF V VAL ( A , VL , UV , A NOl E , M ) 

T M I S SUBROUTINE OETEP M I NES THE FOLLOWING 
VAL Of c OP A: 

A. = I N P 1 1 T VECTOR (SUPPLIED PY TH ,; CALLING PROGRAM) 

VL =MA GNI TI.JDF OF A 

UV=UMT VECTOR PARALLEL TO A 

AMGL c =ANGIF of a measured ccw from x-axts 

D I M F NS I CM Ml),tJVm 
IF (”-100 12, 3 ,4 
I F ( m— 11C) C ,6,7 

COMPUTE VECTOR LENGTH ONLY 

VI =SOPT( A( 1)**?+A(2)**2) 

P F TUP N 

COMPUTE ANGLE OF VECTOR ONLY 

ANGLF = ATAN2( A(2 ) , A ( L I ) 

PFTUPM 

COMPUTE BOTH LENGTH AND ANGLE OF VECTOR 

VL=SOPT(A< 1 ) **2+A ( 2 1**2 ) 

GO TO 2 

GENEPATF UNIT VECTOR ONLY 

VL=EQPT(A( 1 )**2+A(2)**2) 

00 1 1=1.2 
UV ( T ) = A ( i ) /VL 
RETURN 

COMPUTE anglf and generate unit VECTOR 

A NO l E = ATAN2(A(2),A(1) ) 

GO TO ft 
FMO 
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***************************** 

* RANDOM SAMPLING * 

* NORMALIZED CONVOY ROUTE * 
***************************** 



C 



c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 



COMMON/ONE/SGN*YB 



DIMENSION UV< 2) » W1V( 2 ) , W2V( 2 ) . VI V< 2 ) 
1RNV(2I»UW1V(2) ,UW2V(2) » TEMV ( 2 > , WV < 2 1 



. V?V ( 7 I . ft 1 



t 



REAL NORM 

INITIALIZE URN GENERATOR 
POSI T=UP N( 0 ) 

PI =3* 141 5926536 
SUB SPEED 



V=10. 

SAMPLE SIZE 
JKL*1000 

NO, OF ANGLE ( ALFA) INCREMENTS 
IAI=17 

NO. OF RHO INCREMENTS 



I BI=9 

DEFINE BOTH SUB VECTORS 

CALL VINS( V1V,0. »VI 
CALL VINS(V2V,0. t-VI 

SIZE OF DETECTION RAOIUS (DEFINITE RANGE LAW! 
R=30. 

LENGTH OF SUB PATROL 
PL=1 20, 

INITIAL X-POSITION OF PATROL LINE 
XP=R 

WRITE(6,600IR,PL,XP,V»JKL 
INITIAL ANGLE 
ALF=0.0 

LOOP TO INCREMENT ALFA 
DO 800 IA*1»IAI 

INITIAL SHIP SPEEO 
U=10. 

LOOP TO INCREMENT RHO 
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pp on i p= 1 » i b i 



C INITIALIZE LEG COUNT 

si rr=r. 

f rc A= rn<: ( /'L p ) 

SINA = HN( At F) 

C OFF INF VECTOR OF SHIP VELOCITY 

CCLl VT*'S< UV,U*COSA f U*SINA) 

COMHUTF SHIP/SUR SPEED PATIO 
«?HP=U/V 

INITIALIZE SUB DIRECTION FOP DETERMINISTIC 
C Al CULA T TON 

c r,M = l. 

INITIALIZE DETECTION COUNT 
PL T=c. 

CPyPUTF MAXIMUM TIME for calculations 
T m a, X = ( X P + P ) / (IJ*C0SA) + 1. 

PFTEPMINE RELATIVE SPEED VECTORS AND VALUES OE 
MAGNITUDES AND ANGLES 

CALL VADP( V1V,UV, W1V,-1) 

CALI VADD( V2V,UV,W2V t -l ) 

CALL WAL(W1V,W1,UW1V,BW1 till) 

CALt VVAL( W2V,W2,UW2V t BW2,lll ) 

normalize convoy track eor EACH UNIFORMLY 
DISTRIBUTED SUBMARINE STARTING POSITION 
SAMPLE SIZE FOR NORMAL IS DETERMINED BY THIS LOOP 

PP 7 JJ»1.1000 
YP=NOPM( 2C. )-( PL/2. ) 

DETFPMINE IE DETECTION OCCURS 

EACH COUNT OF THE LOOP INDEX REPRESFNTS ONE SAMPLE 
DP fl J=1,JKL 

PETEPMIMF sub POSIT (UNIFORMLY DISTRIBUTED) 

INPUTS 

XPt X-DI STANCE TO PATROL LINE 
YF: Y-COORDINATE OF LOWER LIMIT CF PATROL 

LINE 

PL: LENGTH OF PATROL LINF 

OUTPUTS 

pr: RANGE TO SUB * S INITIAL POSIT 

BO: BEARING OF .SUB • S INITIAL POSIT 

YB (VIA COMMON/ONE/): INITIAL Y-C00ROINATF 

SUB 

P ! V : INITIAL SUB POSITION VECTOR 
CALL WYFP( XP,YP,PL,RO,BO t RlVI 
r TIME TO END OF FIRST (RANDOM) LEG 

IF(SGN)7fC,l05,702 

C T I M r TO GO FCR SUB TO REACH LIMIT OF PATROL LINE 

C IF T N I T Y A I. SUB VELOCITY IS "DOWN" 
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'’ r 0 T ] = ( YP-YF ) /V 
CP TP 70? 

T I Mr Tf GO FPR SUB TC REACH UPPER LIMIT OF PATROL 
I IH C 

702 Tl=( YP+PL-YR) /V 
-"3 T=Ti 

IHITIA L ! 7 F TOTAL TImf 
T H = C. 

INDICATE THAT the FIRST (RANDOM) L FG IS TO BE 
F XAMI NF P 

MM- ] 

rpMrijTf LENGTH of RANDOM lfg 

or to i c* 

l PGP TP INCREMENT leg COUNT 
10 3 DT 1 r T m= 1 ,mm 

CHICK Tp DETERMINE INITIAL DIRECTION OF SUB 
IF<S0N)KltlCP.102 

TNITI/I. DIP.FCTION IS IJP: USE PARAMETERS FDR THE 
UP I.FG 

101 TW=TW? 

C/LL CPPI ( UWV » UV*2 V ) 

w = w? 

CALI C °P I ( R V » W2 V ) 

CALI. C nF I ( VV » V2 V ) 

GO TO \rf- 

INI T IAL DIRECTION IS DOWN; USE PARAMETERS FOR THF 
DPWN LFG 

173 CALL C n PI (UHV, IJW1V) 

JW= T W1 
W = W1 

CALL CPPKWVtWIV) 

CALL C° c nWrVlV) 

C REVOPSE DIRECTION OF NEXT LEG 

IDS SGN=-SGN 

C COMPUTE ^FLA-TIVE LENGTH O c THIS LEG 

GAIL VM|)t T(UWV, T FMV t TW) 

C Dt T r P M I HE VFCTCR FOR FINAL POSITION ON THIS L C G 

CALt VADD( P 1 V, tfmv.PNV , 1 ) 

CALL C P ACHK ( R l V » R NV » UWV » R M IN , IN f J ) 

C INCREMENT LFG COUNT 

SI FG=SLFG+ I . 

C or TP MI HE I F DETECTION OCCURRED ON THIS LEG 

I F(PKTN-R >3,3,11 
11 TH=TH+ t 
A I F ( T N- T MAX)6»B »° 
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ft IF(MM-mC5,9,10 
1 0 CrNTINUF 

or to 8 

9 T= PL / V 
mm- 1 cocr 
10A TW1=T«W1 
TW? = T *W2 
GO TP ic? 

DCT=r=T+i. 

° CrNTTMUF 
7 C r F' T I Ml IP 

XJKl =JKl *K00 
P=PFT /XJKL 
XL FG=SLFG/X JKL 

WP!TF(ft t fcCl J XLEG,U*RHC,Wl ,W2, ALF,P 
’31 U=lJ+10. 

WFITF(t f FC21 
S02 F r P M A T ( / 1 
TCP MF = ALF + PI/36. 

Sin F PR M A T ( 1 H 1 / // // TISt ’PROBABILITY OF RANDOM ENCOUNTER ', / 
1//T6, 'OF -INITE RANGE LAW* /TB, •DETECTION RANGE: *,F5.l, 
?' MILEF '//T6 f ' LENGTH CE SUB °ATPOL LINE: '.F5.1,' MIL' 

?, 'ES'//TF, 'X-DISTANCE TO PATROL L I NE: • » F5. 1 . ' MILES'// 
CT6, 'SUBMARINE S PE EO : • . F5 . 1 . • KNOTS •/ /T6 ,' NO. O c SAM PL ' 
f » ' F S TAKEN POP EACH SCENARIO: ' f l4///T3,' AVG. NO. SH • 

f,'IP SPELT • 

7 f'PFLATIVF SPEFD ALFA PROB. OETECTION'/T5, 'LEGS • 
p,' S PE C D RATIO UPLEG DCWNL EG ' / / 1 
nOl F0RMAT(T4,F5.2 ,?X,3( 1X,F5.1 » , 2X , F5 . I , 3X , F 6 . 4 , 7X, F6. k ) 
103 STOP 
END 



SUBROUTINE WYEBI XP,YP,PL,RO,BO,A) 

I NPUTS 

YP=Y-POSI TION OF SUB ON PATROL LINF 
YB IS UNIFORMLY DISTRIBUTED 

SGN=INIT!AL DIRECTION (EQUALLY LIKELY) OF SUB 
VELOCITY 

CC MM ON /0N C /SGN » YB 
DIMENSION A( 2 ) ,B(2» 

C DETERMINE SUB POSIT RANDOMLY 

YB = YP + I IP N ( 1 ) *PL 

P DETFRMINE INITIAL DIRECTION 

i IF(UPN(n-P.5)5,6,7 
F GO TO ^ 

C UP 

3 SGN= 1 . 

GO TO « 

C DOWN 

7 5GN— 1. 

P DETFPMINE INITIAL RANGF AND BEARING TO SUB 

B A ( 1 ) =XP 
A ( 2 ) = YB 

CALL WAL( A,R0,B t B0,10l) 
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RFTURN 

END 



C 



C 

C 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



SUBROUTINE CPACHK ( R I , RN,UWV,RCP A , IN , J ) 

DIMENSION Rl(l ) ,RNC1) ,UWVU) , X< 2 ) , Y < 2 1 ,CPA < 2 ) , Z ( 2 I 
COMMON/TWO/Y ,N 

R1=REFERENCE VECTOR 
RN=TERMINAL VECTOR 
UWV=UNI T VECTOR OF DRM 



COMPUTE DOT PRODUCT: R1*UWV=X 

CALL VADD(R1»UWV,X,0) 

IF(X(im,2,2 

CPA WILL OCCUR ON THIS LEG: COMPUTE THE DOT 
PRODUCT: RN*UWV=X 

1 CALL VADD(RNfUWV»X t O) 

IF(X(1) )3,3,4 

CPA WILL NOT OCCUR ON THIS LEG: CPA=R1 COMPUTE 
CPA 

2 CALL VVAL(F1 , RCP A , X t D ♦ 1 00 ) 

RETURN 

CPA WAS NOT REACHED ON THIS LEG: CPA=RN: COMPUTE 

CPA 

3 CALL VVAL(RN,RCPA,X,D,100) 

RFTURN 

DETERMINE CPA: ASSUME LEG HAS LENGTH< 250 MILES 

4 DO 5 1=1,500 

USE HALF-MILE INCREMENTS 
P=I/2. 



INCREMENT DRM BY 1/2 MILE 
CALL VMULT ( UWV » Y , P ) 

ADD TO REFERENCE VECTOR 
CALL VADDCRl ,Y,X,1I 

COMPUTE DOT PRODUCT X OOT UWV=Z 
CALL VADD (X,UWV,Z,OI 

IF Z CHANGES SIGN, OR BECOMES ZERO VECTOR X=CPA: 
OTHERWISE TRY AGAIN 



iF(ztm 

CONTINUE 



5,6,6 

DRM LEG>250 MILES 
WP ITF ( 6 ,60 ) I N, J 

60 F0RMAT(//10X, ’LEG' ,1X,I4,1X, , 0F SAMPLE ' ,I5,1X,'IS 
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1GFEATEP THAN 250 MILES'//) 
RF T UPN 

5 CALL VVAL ( X t RCPA,Y,0, 100) 
PFTURN 
EM) 



Sl'PR ri)T I NF VMULT(A,B,X) 

C THIS SUBROUTINE MULTIPLIES THE VECTCP A BY THE 

C SCALAR X TO FORM THE VECTOR B 

0 I ME NS I ON A ( l ) • 0 ( 1 ) 

DO 1 1=1,2 
1 BI I I = X*A ( I ) 

P F T UP N 
END 



SUBROUTINE VVAL ( A , VL ♦ UV, A NGLE , M ) 

THIS SUBROUTINE DETERMINES THE FOLLOWING 
VALUES PE A: 

A= I NPUT VECTOR (SUPPLIED BY THE CALLING PROGRAM) 

VL=MAGNI TUDE OF A 

UV*UNI T VECTOR PARALLEL TO A 

ANGLE*ANGLE OF A MEASURED CCW FROM X-AXIS 



DIMENSION A(l) ,UV(1) 

IF (M- 100 ) 2 ♦ 3 ,4 
A IF(M-1 1015,6,7 

COMPUTE VECTCR LENGTH ONLY 

3 VL=SORT(A( 1)**2+A(2)**2) 

RETURN 

COMPUTE ANGLE OF VECTOR ONLY 

2 ANGLE =ATAN2( A ( 2 ) ♦ A ( 1 ) ) 

RETURN 

COMPUTE BOTH LENGTH AND ANGLE OF VECTOR 

5 Vl=$QRT(A(l)**2+A(2)**2) 

GO TO 2 

GENERATE UNIT VECTOR ONLY 

6 VL=SQRT(A(1)**2+A<2)**2) 

DO 1 1*1,2 

1 UV(I )=A(I )/VL 
RETURN 

COMPUTE ANGLE AND GENERATE UNIT VECTOR 

7 ANGLE=ATAN2< A ( 2 ) ,A( 1) ) 

GO TO 6 
END 
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SUBROUTINE VIN$(A,X,Y) 



THIS SUBROUTINE INSERTS THE X AND Y VECTOR 
COMPONENTS INTO ARRAY A 



END 



SUBROUTINE COPI(A,B) 

THIS SUBROUTINE GENERATES THE VECTOR A SUCH THAT 
FOR A GIVEN VECTOR B, A=B 



SUBROUTINE VAOD ( A , B ,R ♦ N) 

THIS SUBROUTINE DOES THE FOLLOWING VECTOR 
OPERATIONS : 

IF N>0 : R=A+B 
IF N<0: R=A-B 

IF N=0 : THE 00T PRODUCT OF A AND B IS FORMED AND 
RESULT IS PLACED IN THE FIRST MEMBER OF ARRAY R 

DIMENSION A(l)tB(l)*R(l) 

IF(N)1,2,3 

C ADD VECTORS 

3 DO A 1=1.2 

4 RU)=A(I)+B(n 



2 DO 6 1*1.2 
6 Rm=A(I)*B(I) 
R(1»=R( 1»+P(2) 
PE T URN 
END 



DIMENSION A(l) 

a(h=x 




l 

RETURN 

END 




C 



L UU D 1= 

5 RU)«AII)-Bm 
RETURN 

COMPUTE DOT PRODUCT OF VECTORS 



RETURN 

SUBTRACT VECTORS 
1 DO 5 1*1,2 



C 
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r. THI'') SUBROUTINE GENERATES N(0 , sigma SQUARED) 

<- RANDOM VARIABLES BY THF BOX-MUELLFP METHOD 

r TWO I NPF°FNPENT NORMALLY DISTRIBUTED RANDOM 

r V A P I A B l c S APE PRncUCrD EACH T J M r THIS SUBROUTINE 

r T S C A l L f D 

r P=S T ALCARD CSVIATIQM 

PS M FUNCTION NCRM(n) 

DATA N,TKPFI/i,A.233iP53/ 

GC TP ( 1.2) ,N 

1 o = SPPT(-2. *ALOG(URNU ) ) > 

ARG = TViPP I * ( URN ( 1 ) ) 

NrPM=o*C PS ( ARG ) *0 
X?=P*SIM(APG)*C 

N = ? 

Rf TURN 

2 Mp.RM = x? 

N=1 

cr t^pm 

RNP 
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APPENDIX E 



A. GEOMETRIC INTEGRATION 

Figure 7 is a graph of a family of curves plotted from 
Table E-la which represents the conditional probability 
that detection occurs as a function of the crossing point 
displacement from the patrol line center. The normal 
density function, with mean zero and standard deviation 20 
miles, is also plotted in the figure. The upper set of 
numbers along the abscissa represents the normalized values, 
i . e . , 

x' = x/20 

where x' is the normalized value of x. The product of the 
normal density curve with each of the conditional proba- 
bility curves is tabulated in Table E-lb, and plotted as 
a family of curves in Figure 9. If a grid of 0.1 milli- 
meter squares is superimposed on this figure, the area 
represents one-half of the integral of the function. The 
area under each curve was determined by the grid-counting 
method, and was also determined by use of the Amsler Polar 
Planimeter. Each curve was traced five times with the 
planimeter stylus, and an average planimeter dial reading 
taken. The results of these two methods of geometric 
integration are tabulated in Table 9. The average plani- 
meter dial readings were normalized to 0.998 (the value of 
P(D) for p = 0) using the following equation: 
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TABLE E-l 



CONDITIONAL PROBABILITY OF DETECTION 
GIVEN CROSSING POINT IS DISPLACED 
Y MILES FROM PATROL CENTER 



a. NO ASSUMPTION ON DISTRIBUTION OF CROSSING POINTS 



y 




P(D|Y p = y) 

SHIP- SUBMARINE oPEED RATIO (p) 








o 

II 

Q_ 


P-1 


P-1.5 


P»2 


P = 3 


p=5 


p = 9 


P=°o 


0 


1 


1 


.9354 


. 8000 


.6250 


. 5104 


.5021 


.5000 


5 


1 


1 


.9260 


.7844 


.6250 


.5260 


.5021 


.5000 


10 


1 


1 


. 8687 


. 7635 


.6177 


. 5469 


.5021 


.5000 


15 


1 


.9781 


. 8062 


.7125 


.6104 


.5604 


.5292 


.5000 


20 


1 


.9365 


.7437 


.6500 


.5969 


.5604 


.5292 


.5000 


25 


1 


.8844 


.6812 


.5875 


.5552 


.5469 


.5292 


.5000 


30 


1 


.8219 


.6187 


.5302 


.5135 


.5052 


.5010 


.5000 


50 


1 


.5719 


. 3844 


. 3635 


. 3469 


. 3385 


. 3344 


. 3333 


70 


1 


. 3219 


.2177 


.1969 


.1802 


.1719 


. 1677 


.1667 



b . 


NORMAL DISTRIBUTION (STANDARD DEVIATION 20 
OF CROSSING POINTS ASSUMED 


MILES) 




y 


<t> (y) * 


SHIP 


p ( d 1 Yp = y) 

SUBMARINE SPEED RATIO ( p) 








O 

II 

CL 


P-1 


p=l . 5 


P-2 


P=3 


P=5 


p=9 


P = oo 


0 


.3989 


. 3989 


. 3731 


. 3191 


.2493 


.2036 


.2003 


.1994 


5 


. 3867 


.3867 


. 3581 


. 3033 


.2417 


.2034 


.1942 


.1934 


10 


. 3521 


. 3521 


. 3059 


. 2688 


.2175 


.1926 


.1768 


.1760 


15 


. 3011 


.2945 


.2428 


.2145 


.1838 


.1687 


.1574 


.1506 


20 


.2420 


.2266 


.1800 


.1573 


.1444 


.1356 


.1281 


.1210 


25 


.1826 


.1615 


.1244 


.1073 


.1014 


.0999 


.0966 


.0913 


30 


. 1295 


.1064 


.0812 


.0687 


.0665 


.0654 


.0649 


.0648 


50 


.0175 


.0100 


. 0067 


.0064 


.0061 


.0059 


.0058 


.0058 


70 


.0009 


.0006 


.0002 


.0002 


.0004 


.0002 


.0002 


.0003 



*Extracted from CRC Standard Mathematical Tables, 16th ed., 
1968, pp. 549-555. 
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Area under Curve 



(0.9980 x average reading)/762 (2) 



where the constant 762 is the average planimeter dial 
reading obtained by tracing the curve for p = 0. 

The values contained in Table E-2 for the column under the 
heading "Planimeter P (D) " were computed by equation (2) 
using a Friden electronic calculator. 



TABLE E-2 

RESULTS OF TWO GEOMETRICAL 
INTEGRATION METHODS 



SPEED RATIO 
(( P ). 


P (DETECTION) 




GRID COUNTING 


AMSLER POLAR 
PLANIMETER * 


0 


0.9980 


0.9980 


1 


0.9415 


0.9319 


1.5 


0.7895 


0.7922 


2 


0.6816 


0.6828 


3 


0.5781 


0.5797 


5 


0.5233 


0.5234 


9 


0.5008 


0.5011 


OO 


0.4885 


0.4883 



*The planimeter dial readings were normalized so 
762 increments on the dial represented the known area 



that 

0.9980. 
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n n n n n 



APPENDIX F 



LIMITING CASE FOR P (DETECTION) 



LIMIT 



J=0 
R=30 . 

PL=120 . 

YP=40. 

WRITE (6 ,600) R, PL 
Z=YP+PL 
DO 1 1=1,41 
IF(YP-R) 3,2,2 

2 P=0 

GO TO 4 

3 IF (YP+R) 5 ,5 ,6 

5 P=(2.*R)/PL 
GO TO 4 

6 Y=AMIN1 (ABS (YP) ,ABS (Z) ) 

S= ( YP* Z ) 

W=S/ABS (S) 

P=(R-W*Y)/PL 

4 J=J+1 
IF(J-6) 7,8,9 

8 J=0 

WRITE (6,602) 

7 WRITE (6,601) YP,P 
1 YP=YP-2 . 5 

9 STOP 

600 FORMAT (1H1, 15 (/) ,20X, 'LIMITING CASE FOR P (DETECTION) '/ 

1//5X, 'DETECTION RADIUS (DEFINITE RANGE LAW): ',F4.0//5 

2X , 'PATROL LENGTH: ' ,P4 . 0///7X, ' YP ' , 8X, 'P (D) '//) 

601 FORMAT (5X,F5.1,5X,F6.4) 

602 FORMAT (/) 

END 
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TABLE F-l 



LIMITING CASE FOR P (DETECTION) 

DETECTION RADIUS (DEFINITE RANGE LAW): 30. 
PATROL LENGTH: 120. 



YP 


P(D) 


40.0 


0.0 


37.5 


0.0 


35.0 


0.0 


32.5 


0.0 


30.0 


0.0 


27.5 


0.0208 


25.0 


0.0417 


22.5 


0.0625 


20.0 


0.0833 


17.5 


0.1042 


15.0 


0.1250 


12.5 


0.1458 


10.0 


0.1667 


7.5 


0.1875 


5.0 


0.2083 


2.5 


0.2292 


0.0 


0.2500 


-2.5 


0.2708 


-5.0 


0.2917 


-7.5 


0.3125 


-10.0 


0.3333 


-12.5 


0.3542 


-15.0 


0.3750 


-17.5 


0.3958 


-20.0 


0.4167 


-22.5 


0.4375 


-25.0 


0.4583 


-27.5 


0.4792 


-30.0 


0.5000 


-32.5 


0.5000 


-35.0 


0.5000 


-37.5 


0.5000 


-40.0 


0.5000 


-42.5 


0.5000 


-45.0 


0.5000 


-47.5 


0.5000 


-50.0 


0.5000 


-52.5 


0.5000 


-55.0 


0.5000 


-57.5 


0.5000 


-60.0 


0.5000 
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1.0 



0.9- 

0 . 8 - 



0.7- 



Ship-Submarine Speed Ratio: 00 

Submarine Patrol Line Length: 120 

miles 

Detection Range (Definite Range 
Law) : 30 miles 



0.6 * 
P (D) 




Crossing Position in Miles (from Patrol Line Center) 



Figure F-l. Probability of Detection vs. Crossing 
Position (Perpendicular Tracks) for the 
Limiting Case 
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Naval Operations, Study 548, The Effects of Increased 
Merchant Ship Speed on Vulnerability to Submarine 
Attacks (U) , p. 31, February 28, 1955. (C) 

2. Operations Evaluation Group, Office of the Chief of 

Naval Operations, Summary Report 10, Countering the 
Threat of the Nuclear-Powered Submarine (U) , p. 3, 
April 30, 1962. (S) 

3. Center for Naval Analyses, Operations Research Group, 

Report ORG(M) 44, An Analysis of U.S. Coastal Trade 
Convoys (U) , p. 7 , December 18 , 1943 . Tc) 

4. Operations Evaluation Group, Office of the Chief of 

Naval Operations, Report No. 56, Search and Screening , 
by B. O. Koopman, p. 106-107, 1946. 

5. Ibid , p. 107, Figure 8. 

6. Beckenbach , E. F. , Modern Mathematics for the Engineer, 

p. 279, McGraw-Hill, 1956. 

7. Operations Evaluation Group, Office of the Chief of 

Naval Operations, Study 647, The Vulnerability of 
Independent Merchant Ships (U) , by Neufer, J. E. , 
Appendix A, June 7, 1961. (S) 

8. Koopman, op. cit . p. 107. 
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